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SUMMARY 


The present paper is one of a series of papers to ex- 
tend the analytical methods used so successfully in the 
theory of an incompre s siUla fluid to the case of a compress- 
ihle fluid. 

A stream function of an irrotational flow of an incom- 
pressible fluid satisfies the Laplace eo^uation and, con- 
versely, the, imaginary part of an arbitrary analytic func- 
tion of a complex variable can be considered as the stream 
function of a possible flow; many results of the highly de- 
veloped theory of analytic functions can thus be interpret- 
ed as theorems on the motion of an incompressible fluid. 

By using the hodograph method (introduced in the theory 
of compressible fluids by Chaplygin) a formula had previous- 
ly been obtained for the stream function of a possible sub- 
sonic compressible flow in terms of an arbitrary analytic 
function of a comulex variable; procedures for using some 
methods and technic^ues of the theory of analytic functions 
in the theory of subsonic flows have likewise been indicated 
and, as a conseo_uence , new flow patterns have been obtained. 
These flow patterns include examples of flows around symmet- 
ric and ncnsymmetric obstacles, under the assumption of the 
true pressure-density relation, p=<?p*^, o and k being 
constant s . 

In this paper the foregoing results are improved and 
completed. A formula (analogous to that for subsonic flows) 
is derived, which represents a stream function of a possible 
supersonic flow in terms of two arbitrary differentiable 
functions of one real variable. Finally, sonie instances are 
discussed in which flow pattern defined in two neighboring 
parts of the plane can be combined into one flow pattern de- 
fined in the combined domain. This last method leads, in 
some instances, to partially supersonic flows. 
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INTRODUOTION 


The development of. research in compr essilile fluid the~ 
ory has made it desirable to have adequate mathematical 
tools for dealing v/ith problems of compressible flows. One 
of the reasons for the success of mathematical methods in 
the study of two-dimensional irrotational steady flows of 
an incompressible fluid is based on the fact that it is 
possible to represent the stream function of such a flow as 
the imaginary part of an analytic function of a complex 
variable. As a consequence, various results in the highly 
developed theory of analytic functions can be applied to 
yield solution^ of problems- in hydrodynamics. 

In previous publications of the author (references 1 
through 5), this approach has been generalized to include 
the case of compressible fluid flows; this was accomplished 
by representing the stream function of a possible flow of a 
compressible fluid in terms of two arbitrary functions of 
one variable. In that part of the flow where the character 
of the flow is subsonic, one of these functions is an ana- 
lytic function of a complex variable, the other, its conju- 
gate. In the region of the flow in which its character is 
that of a supersonic flow, each of these two functions is a 
different function of ono real variable. 

The development of this approach raises several complex 
and fairly difficult questions; the purpose of this report 
will be to discuss, in some detail, the problems entailed by 
these methods . 

In order to facilitate reading, in sections I and II 
the general idea of this method of approach will be summa- 
rized. 

This investigation, carried out at Brown University, 
was sponsored by and conducted with the financial assistance 
of the National Advisory Committee for Aeronautics. 

The author would like to express his sincere apprecia- 
tion for the assistance and advice he received from Mr. 
Leonard Greenstone and to thank Mr. Herman Chernoff for his 
help with section III and Mr. Bernard Epstein for his help 
with section IV. 
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SOTA!TIO!r 


In dealing with differential equations the following 
notation is often used; 


&u 

u = r — = 
z bz 

1 /Bu 

2 

Bu \ 

u- = ^ = i ( 
z Bz 2 ^ 

'p- + i 

<ox 

BuN 
By / 

B^u 1 

+ 


u; z = X + iy. 


- iy 

B z BF 4 

VBx® 

ByV 4 "" 

Z ^ X 

a = [ao^ - 

■1''= 

)>/ a 

- Dq^ 

, speed of sound 



8.Q speed of 

sound at a 

stagnation point 




aj^ coefficients in the series expansion of T in powers 

of also used in the sense of equation (72)” 

hj^ coefficients in the series expansion of in 

powers of ; also used in the sense of equa- 

tion ( 73 ) 

Cjn (See equations (146), (148).) 

d-n ^ ^ ^ 1 An^^\ ^n^^\ (See equations (86), (87).) 

e "base of Naperian logarithms 

exp (x) = e^ 

f, g arbitrary analytic functions, in the subsonic case; 

arbitrary twice differentiable functions of one 
real variable in the supersonic case 

_ dg(z) dg^°^ (z) 

dZ dZ 

g[o] (Z) = g(z) 
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g ] (Z) nth iterated integral of g(Z) 


h 

k 


( 


k - 


k + 



h 


~ for k 

76 


1.4 


constant in the equation of state p = ap^; the 
ratio of specific heats at constant pressure to 
constant volume; k = 1,4 for air 



^ (1 - M®(H)) 


p pressure 

Pq pressure at stagnation point 

q magnitude of the velocity vector 

q = qe ° velocity vector 


schlicht= univalent 
s = 1 - T 

Sn (See equation (87),) 

u, V Cartesian components of the velocity vector 

X, y Cartesian components in the physical plane 

A constant in the Von Karm^n-Tsien equation of state 

p = A + — 

P 

B = for M > 1 


(1 ) /A 

E ' (A) = - / ( Ai) d Ai 

° A 

g(n+i)(A) = ^ r (By^^^^^(Ai) +n(Aj) E^”^(Ai))dAi 


A 
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(See equations (159).) 

£^^)(A) = E^(A*) 

F = - (N. + N®‘) = - ^ ■ [- (Kk - 1) - 4 (3 - 2k) M® + 16 ] 

^ 64(1- M^)® 

for k=1.4 (See equation (l9).) 


F, = Ni‘ 


a .. I SNi 

‘ 2 a ^ 


(See equation (45).) 


Fjj^ polynomial approximation of the mth degree in e* 

to F 


Fg . (See equation (l58).) 


H = exn 


1 

/ 1 1 2(k-l) 

) / Nd(t+t)j=* ^ i 

(l-M^s)^^*^ l+^(k-l)M® 

\_ 45 J 

in this sense H is used only in the series expansion 
of '4/ (See equation (17).) 


dq: in this sense H is used only as an independ- 

ent variable (See equation (32), ) 


= exp^- / lli(Ai)dAi^ 


Im imaginary part of 

D ® 

L(H) = (M® - 1) — (See 1(H).) 

p(H)® 

L*'’^^(2\) =» — H(2^)^^^^ (2A) 

2 nj 
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M local Mach nuraher; M = q/a 

Mj , Ms, M 3 , M 4 (See equation (57).) 


N » - 


Ni = - 




(k+ 1 ) 


M 


S (1-M^)""" 


(k■^ 1 ) 


M 


8 (M® - 1 ) 


3 /a 


q " (X) = -4 / P dX 
•.y —00 

(n + i ) 


Q, 


(X) = - 


2 n + 1 


(4f (^i)q(“)(X 3^) + 'iXX^'^^Xi) )dXi 


(See equation (l07), and appendix II,) 


0,^^^ computed employing Pjjj instead of P 
Re real part of 

S (») = i_C'i Sjtv 0 

P® &0® pSqVp Sqy 

T =/l - M® 

U, Ui, Ug, U*** (See equations (84), (85 ), ff., (173), (174).) 

P(A, 6) = Hj^“^\l/(A, 0) (See equation (50).) V = V 1 +V 3 

00 

ViU, n) = f(|)+ y b^^\a) f (e) 

n=i 
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00 


Vs(t, r,) = g(Ti)+ y Efn)(A) gC^3(n) 


n = 1 

V=^, V 2 *, Va*i‘ (See footnote, p. 24.) 
N(|, n) (See equation (l58).) 

X (See equation (lOO), ff . ) 

Z SI A + i 9 

8 = l/k-l; 8 = 5/2 for air 


j = (-3(1= + 1) I V"" 

■n=A-6 ; ri* = n- a/2 (See equation (l55). ) 

6 the angle the Telocit 7 vector ~Hi makes with some 

fixed direction; "q = qe^®; also 0 - •^( I ~ ) 

(See equation (46).) 


A 


1 

2 


log 


1-7 

1+7 


+ 


1 

h 


log 


h~^ + 7 
h~^ - 7 


(See equation (20).) 


I = A + 0 

I* = I - a/2 (See equation (l55), ) 
p density 

Pjj density at a stagnation point; equation (ll) 

0 = T - 0.15, equation (74); also a constant in the 
pressure-density relation p = ap^ 


T 


T 


(k - 1) 

2ao® 

1 3 


for air, assuming a. 


1 


v{/ stream function 
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\|/* =HiV* 
A 


, & ®cp 6 ®q3 . 

Laplace operator: Atp = — — -j- — — = 4 


(— - ) 

' ft ^ ft 9: y 


H 




5y‘ 


■5z & z 


A(H) = 


Vl(H) dH = ^ tan“^(hB) - tan ^ 

ll 


B 


(V) = + n(g + n) V = 0 

n( 0 an analytic function of a complex variable satisfying 
certain conditions 


Remark ; Observe that auite frequently functions will he con- 
sidered in different planes although the notation will not, 
in general, indicate this. Thus, given f(x, y), let 


J 0 

5.(x, y) 

to obtain f(x(x^, x®), y(x^, x®)) = f^(x^, x®). The super- 
scrint will, in general, be omitted and only f(xi, y® ) 
written, since the meaning will be clear from the context. 


X = x( X^ , X® ) 
y = y (x^ , X® ) , 


I 


In the following, the fluid flow to be considered will 
be supposed to be that of a two-dimensional irrotational 
steady flow of an inviscid, compressible fluid. 

The assumption of the law of the conservation of matter 
leads to the equation of continuity: 

^(pu) + -~(pv) = 0 (1) 

ox dy 

where p is the density, (u,v) the Cartesian components of 
the Telocity vector, and (x,y) the Cartesian coordinates in 
the plane of flow. 

This equation implies that the differential expression 
pudy - pvdx is a complete differential - that is, that it 
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is the differential of some function, say \i/(^iy)i the so- 
called "stream function of the flowc" 

Thus 

dvj/ = pudy - pvdz (S) 

which clearly implies the relations 


5\j/ 

hx 





( 3 ) 


•» 

The assumption that the flow is irrotational may he ex- 
pressed mathematically hy the equation 


du _ 
dy 6x 


(4) 


This implies that udx + vdy is the differential of some 
function, say, cp , the so-called "potential function of 
the flow," Thus, 


dcp = u dx + T dy 


( 5 ) 


« 


so that 


^ = u. ^ = V 
dx Sy 


( 6 ) 


and 


q, the magnitude 

q 

Equations (3) and 

bcp 

dx 


of the Telocity vector, 



( 5 ) imply the system of 
p dy’ dy p 5x 


is given hy 

(7) 

equations ; 

( 8 ) 


relating the stream function and potential function of a 
flow described ahove. 
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If it “be assumed that the fluid motion, represents an 
adiabatic process, then the equation of state assumes the 
form! 


p = a + A ( 9 ) 

where A, a and k are constants (k c l.U for air), and 
p is. the pressure. 

As the equation of state expresses the pressure as a 
function of the density alone, the first integral of 
Newton's law of motion for a gas may he easily obtained in 
the form of Bernoulli's equation: 


1 

2 


qa 



dp 

P(p) 


0 


( 10 ) 


The pressure, p, may be eliminated with the aid of some com- 
putation (which is omitted here) by combining (9) and (lO) 
to obtain 


P 


1 . 1 



k-1 

■ k-1 r /a 

cp Y , 'f 1 ' 

L 3 ao®J 

1 

u 

■D 

0 




(11) 


Here are the density and the speed of sound, re- 

spectively, at a stagnation point. 

If p is eliminated from equation (s) (giving p its 
value in (ll) and substituting for q the expression in 
equation (?)), it may be seen that /8) represents a system 
of two nonl inear partial differential equations for p and 
\ 1 / . 


In the case of an incompressible fluid flow, where p 
is constant, equation (8) represents a system of two linear 
partial differential equations, and it is these equations 
which provide the means whereby hydrodynamics may be studied 
as an application of the theory of analytic functions of a 
complex variable, 

In order, in the case of a compressible fluid, to ob- 
tain a system of linear equations, it is, necessary to replace 
this approach by an alternate and more complex one, namely. 
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the "hodograph method,’^ according, to which cp and are 

■ considered aa functions not of the Cartesian coordinates of. 
the plane of flow, (x,y), hut of q, and 0, where qei£ 
is the velocity vector. ihinction q denotes the speed and 
6 the angle which the velocity vector forms with some fixed 
direction. 

Consider a steady two-dimensional flow pattern; then at 
every point ( 3 r,y) there is a certain velocity vector qe^® . 
Consequently, to every point (x,y) of the flow pattern, it 
is possihle to associate the pair of quantities, q and 0» 
If q and -Q are now considered as the polar coordinates 
of .some new plane, the ''hodograph plane,” to a streamline 
of the given flow pattern, there will correspond a line in 
this plane. The image of .the flow patterns so obtained is 
called the hodograph of the flow pattern. 




In addition 
hodograph plane, 
of the flow in 
the Cartesian 
viously, it is 
tial functions 
planes , 

That is, if (x,y) are 
given as functions of q 
and 0 , 

X = Cq, e) 

y = (q., 6) 

or of log q and 0 


flow in the 
have the image 
is , the plana 
and 0 . Oh- 
and poten- 
logar ithmic 



to having the Image of the 
it is frequently useful to 
the "logarithmic plane,” that 
coordinates of which are log q 
possihle to consider the stream 
as defined in the hodograph or 
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X = X3 (log q., 6) 
y = 7a(lo€ <lf 8) 

then cp, \}/ evidently may he expressed as a function of these 
variables 

cp(^,y) = cp[ 35:1 (q., s) ,yi (q , 6) ] =cp[x3(log q,G),y3(log q,e)] 

\1/ (x,y ) = \l/[ xi (q, 6) ,yi (q, e) 3 = ^[ X3 (log q, 0) ,ya (log q, s) ] 


On the other hand, if the potential and stream functions are 
given in the hodograph or logarithmic planes, that is, if 


and 


cp 


cpi (q, 6 ) , =* 'ka (q, 0 ) 


cp = cp2(log q,e), \ 1 / =r \t^ (log q,Pi) 


respectively, then it is not difficult to determine the 
corresponding streamlines in this physical plane (the actual 
plane of the flow). 

In the ca'se of an incompressible fluid flow, the po- 
tential and stream functions considered in the logarithmic 
plane, that is, as a function of log q and 0, satisfy 
a system of partial differential equations which has exact- 
ly the same form as equations (s). Thus, 


lil * a ^1^ 

50 pd (log q) ^ 


a CP 

a log q 


p ae 


(p being constant). 

Therefore, except in problems of a special character, 
it is more convenient, in the case of an incompressible 
fluid flow, to operate in the physical plane (the (x,y)- 
plane ) . 

The situation changes completely in the case of com- 
pressible flows, for, as has been mentioned, cp . and \1/ 
considered as functions of x and y satisfy a system of 
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nonlinear partial differential eauations. If, however, m 
and ij/ are considered as functions, not of x and y , out 
* of log q and 0, than, as Molenhroek and Chaplygin have . 

shown, *cp and .vK satisfy a system of linear partial differ- 
ential equations : namely. 


dcp _ 

5G p &q 

Sep _ 1 - 

Sq pq S9 


V 

f 




( 12 ) 


where p is given hy (ll). (See reference 6. ) 

If cp is eliminated, then the following partial dif- 
ferential equation of second order is otitained for \1/: 

S-(>t/) a + 3. 0 (IS) 

p® pdq'^pSq'' 


In the case of an Incompres sihle fluid, So(\{/) 
an ap-proptiate choice of units 

— ^ + T = 0 

Se® d(log q) 


hecomes for 


(14) 


and, hence, the general representation of a possihle stream 
function in terms of an arbitrary analytic function of one 
complex variable may be obtained! namely, 

(log q, 6) = Im|g(^^ , 6+ i log q (l5) 


where ^ is an arbitrary analytic function of the complex 
variable ^ • 

In previous publications of -the author , a generaliza- 
tion of (15) for the case of a compressible fluid has been 
provided. The formulas of this generalization are the means 
by which various flow patterns of a compressible fluid may 
be obtained. (See eauations (24), (25), . . . ,) 

Naturally, a procedure for obtaining flow patterns is 
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only the first step in the development of the theory, since 
in moat cases it'^.s not "some" flow which is required, hut, 
rather, the flow around a given profile which must he known 
(as well as the laws which govern the motion of a compress- 
ible fluid). However, it may he noted that, in the subsonic 
case, at least, an approximate solution of the problem may 
he obtained if, for the function g in equation (24), the 
function which represents the complex potential in the log- 
arithmic plane of the corresponding incompressible fluid 
flov; is substituted. 

In the following, it will be convenient to consider the 
subsonic and supersonic cases separately. This will be done 
in sections II and III, respectively. 


II - THE SUBSONIC CASE 


As was indicated in reference 3. the function (which, 

when multiplied by 


where 


H 



(1 - T®) 


-j-B/4 

J 


- (1 



1 /4r 



(16) 


T® = 1 - M® 


(17) 


yields the stream function xj/ ), considered as a function of 
^ and b satisfies the equation! 


where 





0 


Z = A + 10 


(18 ) 


F 


-0.12T® + 0.51 + 

T^ 


0. 63 0,45 

T^ T® 


( 19 ) 


and T and A are connected by the relation 
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Expressions (I 6 ) , (l9), (20) have been evaluated for k = 

1.4. In references S and 5 they are given for an arbitrary 

k. 


If E is identically zero, then (I 8 ) becomes the 
Laplace equation 


V 


1 r ^ -6211 

4 6 X® 60 ® 


( 21 ) 


and 


\{r * = Im [g(Z)3 , Z = X+ i6 


( 22 ) 


is the "general solution" of this equation; here g is an 
arbitrary analytic function of the complex variable Z. 

Remark ; Eor k = -1, E vanishes - that is. in this case 
the "compressibility equation," in the (X,0 )-plane (the so- 
called "pseudo-logarithmic" plane) becomes the Laplace 
equation. 

The Von Karman-Tsien method of obtaining compressible 
flows around closed bodies is based on this assumption (i.e,, 
that the pressure-density relation is of the form 

a 

p = A + — • 

P 

instead of the adiabatic relation (9) employed by the au- 
thor). The stream functions in the (\,0;-plane are then 
harmonic functions. (See references 4 and 5.) 

As was proved in reference 2, the representation (22) 
of the "general solution" in terms of an arbitrary analytic 
function of a complex variable can be extended to the case 
where E is not identically zero, but is any function 
satisfying certain conditions. 

Under the assumption that the function E* and its 
derivatives satisfy the inequalities 
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K 

d X 


c(K + 1) : 

- K+3 ’ 

(-A) 


X < 0 


and E= 0,1,2,. . . (83) 


where c is a fixed constant, it was shown in reference *5 
that the expression 


Im H( 2 X) 


g(Z) + y l||Ii -i(n)( 3 X)gC^l (Z)p ( 34 )^ 


t\= 1 2 n i 


J 


where 


gl^°^(z) = e(z) 




> 


gf^''^(z) = r (zi)dZi 


(35) 


m “ X|3|3|* « • 


(g(Z) an arbitrary analytic function of the complex variable 
Z; represents a solution of 


U _ + F-^U = 0 
ZZ 

0?he difficulty which arises in the further development 
of this approach arises from the fact that the function F 
which actually appears in (I8) is a rather complicated func~ 
tion of A, and it is pot clear offhand whether or not it 
satisfies the inequalities (33). In order to overcome 
this difficulty, two possible alternatives may be employed: 

1. In reference 5 it was proved that in the domain 
A < 0 (which is the domain under consideration in the sub- 
sonic case) F, given by (l9), may be represented in the 
form 


^The 0, depend only on the function "S*' and there- 
fore being independent of the particular function g(Z) can 
be computed and tabulated once and for all. (See appendixes 
I and II.) 
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I = 


CO 


n^Ta 


'n 


2nX 

e 


( 26 ) 


(See ta"ble 2.) 

If, now, the infinite sum (26) is replaced hy the 
finite sum 

. m 

^m = y (27) 


(or approximated hy some other polynomial in e® ) , then, 
as may he shown rather easily (see reference 3), for every 
finite m, Pm and its derivatives satisfy the inequalities 
(23); and hence the representation (24) may he employed to 
yield solutions of (18) (P replaced hy Pm ^ ‘ 

As has been pointed out in previous puhlications of the 
author, the solutions of (l8), when P is replaced hy P^, 
will then assume the following form 


* r 

= Im j^H(2^) 


g(Z) + 


z 


n=a 


(2n)! 

2^’^ni 


<i^(^>(2A)gKl (z) 


1 

J 


(28) 


(n ) 

the suhscriut m in, indicating the dependence of 

(n ) ' “ 

Q' ' on the m chosen, so that it seems liJcely that the 
actual solution of (l8) may he written in the form 




1 i m \1/ 

m 

Bl-> CO 


(29) 


and this was shown to he the case. (See reference 3, ) 

The disadvantage of this approach lies in the fact that 
for values of M close to M = 1 (corresponding to values 
of T near T = 0), it is necessary to use a large number 
of terms of (26) in order to obtain a good approximation for 
P; that is, m of (27) must he chosen rather large. 
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In appendix I of reference 5, the developments of T 
and in powers of ware given. By substituting 

these expressions for T and in (l9), the desired 

expansion of T in powers of e“ was obtained. In this 


appendix the values of that 

CO 

T = ^ ®'ne®’^ , are compared with 

n=o 


la , 


I 

\ = r 


&^e 


an^ 

» 


where 


the exact values of T, 


and the difficulty mentioned above for high Mach numbers may 
then be observed. 

In this case the determination of the Q,jjj ^ ^ does not 

involve any theoretical difficulties, but for large m, 
does entail difficulties of a computational nature. 


Bemark : Note that in the neighborhood of Q? » 0 , that is, 

M = 1, I and l/T can be developed in a series of 
i/s , 

S = (-3.6A) .1 If these series are substituted in C19) 

for T and l/B,- a development of S' in a power series of 
I, which holds for M < 1, is obtained in the neighborhood 
of M «= 1 . 

3. The second alternative consists in proving that the 
function P as defined by (19) does satisfy the inequali- 
ties (S3), and therefore by the theorem stated in appendix 1, 
the series (24) will converge. 

Nquat'lons ( 96 ) determine each up to a constant; 

this constant can be specified by the requirement that foT 
some fixed value, say q = q^ (and, therefore, for some A, 

say, A = ^q), vanishes. If the same q^ is chosen 

for all n, n = 1,2, . . then the physical meaning of 

such a choice is the following! 


Consider a flow, the speed of which at every point is 
nearly q^; then, in the pseudo-logarithmic (the (X.,6)- 

plane), the flow resembles that of a distorted flow of an 
incompressible fluid. - 

^Jor general k, ^ = (~3(k+l) (X/2) ) ’•/^ . This. choice 
yields T ^ . . 

®Note that ‘in the incompressible case, X and log q 
are identical. 
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The most natural choice, then, to make for would 
he^ Qq = 0; that is, assume all vanish for X, = — 

In appendix I, proof is given that the series (sU) con- 
verges, under rather general conditions on F which include 
in particular the case defined hy (I 9 ). 


For values of K close to 1, many terms of (24) have 
to he taken into account in order to obtain good approxi- 
mations to the stream function; oonaequently , it is neces- 
sary to have as many of the coefficients 


(n) _ (2n)l 


O « 


(2X) 


(so) 


of 




= InV y 




X.(n)^[n] 


n =0 


(31) 


In appendix II, expressions for 


(2X), 


as possible. __ v 

n * 1,3, ... ,8 are derived, the first four of which can 

be expressed in closed form as a function of T, while In 
the expression for the last four, an integral appears; how- 
ever, the integrand in this case can be written in closed 
form, so that at least numerical computation of the Integral 
will not prove too difficult. 


Tables of several of the functions involved have also 
been computed. (See table 5.) 


^Except for = 0, there is no reason to distinguish 

any q^ from any other, and if q^ be chosen different 
from zero the question arises as to why this particular value 
0 ^ Iq ws-s chosen and not some other. 
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III - THE SUPERSOHIC CASE 


Every function \j/ which satiefiee equation (13), or 
(setting 


the equation 


9 



0 


■'^ee '!^hh = 0 


(32) 


(33) 


can be considered as the stream function of a possible com- 
pressible fluid flow. 

In this section M will be taken greater than 1, 

Since both M and p are functions of q, and since 
equation (32) may be interpreted as defining q as a func- 
tion of H; consequently, the coefficient of >l/g0 in (33) 
is also a function of H, say, I>(H) - that is, 

L(H) = r ^ 0 (See Notation. ) (?4) 

P® 

The reduction of equation (33) to canonical form may be ac- 
complished by introducing the new variable AxA(H), de- 
fined by 

H 

A(H) = J y'li(H) dH (35) 

Thus 1 

*hh =.A(iT Ia H)= KhJ'I'aa * 


This H is not to be confused with the H defined by 
eouation ( 16 ). The lower limit of integration in (35) is to 
be chosen later. 
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so that the differential equation (S3) assumes the form 




- 


1 

.2 L(H) 


BL(H) ' 

3A - 



0 


the equation is then in the required canonical form. 


(S7) 


The coefficient of \1 /a, while it is a function of A, 
is difficult to express expiicitly as such, hut may, on the 
other hand, he given comparatively simply as a function of 
M. If SiCSA) is defined hy 


the expression 


1 bL(E) 
2L(H) dA 


4 Ni(2A) 


(S8 ) 


Ni(2A) 


k 4- 1 li 


(39) 


may then he obtained hy formal computation; further A(M) 
then may he seen to assume the form 


A(M) = — tan“^(h*/M® - 1) - tan**^ (</m® - l) 
h 


= i tan ^ (h 


B) - tan 


-1 


B 


(40) 


where 


B® = M® - 1 


(41) 


and 


h 



k - 1 
k + 1 


It is uossihle to simplify equation (37) hy the use of the 
following transformation. 
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\l/( A . e) = Hi V( A. 9) 


( 4 .?) 


Hi := 


(- I Ni(Ai)dAi) 


( 43 ) 


V ( A , 0 ) then satisfies the eauation 

^00 - ^AA ^i(2A) V = 0 

where 


( 44 ) 


I-iCbA) = Ni® + i ^ 

2 ^ A 


X.+ 1 


64 L B 


5(k+l) 12k^ 6k>-14 > , / s 

+ — r-+ ( 4k + 8 ) — ( 3k- 1 )B (45) 


B’ B‘ 


finally, it is convenient to use an alternate form for 
equation (44), If | , q are defined "by 

I = A @ , n = A - 0 ( 46 ) 

then ( 44 ) assumes the form 


+ ^i(l + n) V = 0 ( 47 ) 


The general solution of eauation (47), which is of hyper- 
bolic type, may he represented by a formula which involves 
two arbitrary, differentiable functions of one real variabTe. 

In reference 5, appendix HI, two different methods 
were considered of obtaining such a general solution. The 
first method is essentially based on obtaining the repre- 
sentation by use of P.iemann's function of equation (47), 
while the second is the analogue, for hyperbolic eouations, 
of the representations (84) and (3l), for solutions of the 
elliptic eauation (18). (See appendix III of the present 
r eport . ) 

The validity of this representation follows at once 
from the following; 
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Theor em : Let 


Z'^(v) = + n(| + Ti) T = 0 


( 48 ) 


where ( ^) ( Q is considered as continued for complex val- 
ues of the argument) is an analytic function of the complex 
variable ^ , which function is supposed regular for 




then, there exists a set of functions 




1^2^* • » • 


= L 


(49) 


so that 


V ( I , Tj) = ^ s( I » "H) 

is a solution of (48), where 

CO 

Vi(I.T)) = f(0 + ^ E^’^^(A)f ^°-^ ( O 
* 

CO 

. Vsd.Ti) = g(-n) + y i^“^(A)gI^^^('n) 

n = i 


(50) 


(51) 


and 

= 

f Cn+i](l ) 


fU). 



f li)d 


» 


0 


g ^ ° ) = 


g(Tl ) 



^^^ ( n 1 ) dT| 3, 


o 


f and g being two arbitrary, twice differentiable func- 
tions of I and Ti , respectively. 


The proof of this theorem will be given in appendix III 
Note that in slightly different form this theorem was enun- 
ciated without proof in reference 5. 
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This theorem cannot, however, he applied directly to 
equation (^7) as given hy (4f>) has a pole for A = 0. 

There are, however, several t)os e ihilit i es for overcoming 
this difficulty. - For example, F^_ can he approximated in 
the supersonic range, which is under consideration hy a 
function F^^, which satisfies the conditions on 0 in the 
theorem. 

In many instances, the whole flow lies in a part of 
the supersonic region and the theorem then can he applied 
directly, merely hy shifting the origin, as is done in the 
following example. In appendix III, this will he discussed 
in more detail. 

The following extremely simple example is intended to 
serve only as an illustration of the_,method described above 
for obtaining flow patterns in the hodograph (or'related) 
planes. A discussion of the procedure necessary for deter- 
mining the .corresponding flow in the physical, the (x,y), 
plane will then conclude the section. 

Choose for the f ( & ) , g( P ) of equation (5l) the func— 
t ions'i 

f(|*) » 100(|* + O.l)® 
g(p*) “ — 100(p* + O.l) 

If in (4S) ri = 0 , the corresponding flow in the 
, p)--olane will he 

V* = 100 [(4* + 0.1)=' - (n* + 0.1)®] (53) 

These streamlines are indicated in figure 

In the following the general case Q ^ 0 will he 
cons ider ed . 

The evaluation of the functions 7* , 7 * ® introduced 

iSee appendix III. Essentially, i* , p* differ from 
I , P, respectively, only hy a constant and 7* is the 
function V of- equations (48), (50) when the corresponding 
value of I*, T)* has been suha.tituted for p , rospec— 

t iv ely . 

®See appendix III. Functions 7*. 7| are obtained from 
Vg, respectively, of equation (5l) when the corresponding 
values of | , P* have been substituted for t, p, respec- 

tively. 
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in equations (6l) then will consist of the following stepsi 
f(|*) = f^°^(|*) = 100(1* + O.l)® 


f [^3(1*) 

f [2] (I * ) 



U* + 0,1)® 


+ 0,1)* _ (O.lf 1* _ (O.lf 

3 12 


] 


f[3] (I *) = 100 


(I* + 0,1)^ 

60 


( 0 . 1)3 |*2 ( 0 . 1 )“^ I * ( 0 . 1 )' 


12 


60 


f I*) 


= 100 


+ 0.1)“ (0,1)=^ e 


,3 »*3 


,4 **2 


360 


18 


(0.1)" I 

24 

(0.1)5 I* (0.1)' 


60 


360 


] 


d,f* r T 

In table 6, the values of — — and f^-^(t*)i u =» 0,1, 2, 3, 4 

dt 

for a given set of values of |* are tabulated. 

Remark : Since gCp*) = — f('n*), the values of g^’^^(ri*) and 


dp’ 


may also be obtained from table 6. 

The values of T*, Vg, and 7* = V* + Tg as well as 
>j/* z= H V*, are tabulated in table 10; U* is a function 
of i*, p*. 


In order to carry out the transition to the physical 
plane, the values of d\l;*/oq and Svl/*/9S are needed.^ 

^Henceforth the asterisk which indicates the dependence 
of the function on , p* will be omitted, as this will, 
in general, be clear from the context. 
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4c - n*) + EYp* = H 

SA ^ 


(k+l) 

S 


T + Yfe ♦ 


= H 


( k+l) (BS+i)g 
8 B 3 


V + 


] 


( 56 ) 


where 


i* 




n=i 


(n) ,[n](j.) 

^ 2 ” A* 


n=i 


CD 


11=1 


(5 7 ) 


= Mj_ + Mg + M3 


V 4. = 
T| 


g 


3 'n" 




n=i 


+ Mg + Mg 


and 


: He 

I 

a 


AA B t in _ T 

dq q’ 


’iq" = r . ’ig = -* 1 


( 58 ) 


1 in 

A 



^(n) 

^A’^ = ~SA^’ 


the use of S 


is to call 


attention to the shift of the origin; , Mg, Mg are equal 

to the expression in the first, second, and third brackets, 
respectirely. 
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It follows that 




B 




( k+l) (B^+l) 


B3 


V + Ml + 2Mg + -2M_ + M, 


(59 ) 


\|/fl = H(V 


- Y^*) = H [Ml - M4] 


M - 4 . 


E 


(n) ^[n-i] 


g 


(ti*) 


( 60) 


n=i 


In table 9 the ralues of Mi, Mg, M 3 , and \j/H are 

B H ® 

tabulated , 


In addition, in table 9 the values of 


Ui ) “ ® -1 

\l/= const P'1 


, p S 2 

\l/^q - B ^l/ 

rt a 




( 61) 


and 


(g) 


— s in 0 


\|/ = const 


pq‘ 


1 2 3 , a 

Vl/^q - B te 

— 


are also g iv en . 


As has been shown in reference 5, the corr esp ondinc 
values of x and y, that is, the image in the physical 
plane of the hodograph flow, may be obtained from the formu- 
lasj; 

The integrals are to be understood as taken around the 
corresponding contour. 
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1 1 

— \{/- cos © — \{/n s in 6 d.6 +' 


\J/ „ c o s 0 — ij/ s i n 0 do 

L Pq.^ ® pq. q J 


1 1 1 
y= — \!/ sin 6 + — Uf c cos 9 d9 + 

J L p ' pq ' “ J 


( 62) 


^ 1 

•d/ftSin a + \l; cos 0 do 

pqs ® pq q 


( Pq is supposed set eq_ual to 1 here). 
Now, along a streamline 


d\l; = dq +v|/0 d0 = 


(63) 


that is , 


de = - ^ d. 


(64) 


so that if the integration is carried out along a streamline, 
( 62 ) then assumes the form 


3C = / COS e J- 

./ -el- q 


• rtl 22 ■} 

s in 0 ® '4^0 - q 4/ „ — ^ d-q 

• 0 L 


qjpq' 


( 65 ) 


Similarly, if the integration is carried out along Ti* = 
constant, from the fact that 
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I* - = 2B 

I* + T\* = 2A 


it is easily seen that 


« 1 

dQ = - d£ = dA 

^<1 . q 

da = — dA = — d9 

dA B 


> ( 66 ) 


and therefore along ti* * constant 


X 



B in 0 d^( 

pq 




y 




+ — " cos 0 d\l/ 

pq 


J 


( 67) 


(See diagram II, p. 85.’) 


IV - THE MIXED CASE 


To every analytic function, g, of a complex variahle 
Z there corresponds a solution \]/*(Z) = ^’'‘(A,’0) of equa- 
tion (18) obtained by substituting g into the operator 
given in equation (24). Referring to equation (I 6 ), 

\|/ = \|/*H may then be interpreted as a stream function, in 
the psuedo— 1 ogar ithm ic plane, of a possible flow pattern of 
a compressible fluid. ^ 

^It is necessary to speak of "possible” flow patterns, 
since it may happen that the image in the physical plane of 
a flow obtained in this way, is multiply covered, in which 
case the complete flow pattern has no physical significance. 

( It should be noted, however, that, in this case, those parts 
of the flow pattern which e "schlicht"” do possess physical 
s ignif icance . ) 
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This flow pattern is defined in every (simply connected) 
intersection of the domain of regularity of g, with that of 
the domain L (L is that domain of the (X, 9)— plane for 
which 6®< 3X®, \ < O), which int ei-s e c t ion is supposed to in- 
clude the origin. 

If, for g, the complex potential of an in compr es s ih 1 e 
fluid flow (in the pseudo— logarithmic plane) is substituted, 
the resulting function is defined in a domain . H which in 
many instances may lie partially outside L. (See fig. 5.) 

e 


X 


S'igur© 5. 

Despite the fact that the function , may be regular 
in the whole domain H, the procedure which has been de- 
veloped so far, yields only those values which lie in that 
part of S which lies in L, so that the auestion arises 
as to how to decide in what regions outside of L, xj/ * ex- 
ists, as well as the question of how it may be determined. 

In order to evaluate for those values of X which 

are near zero — that is, for high Mach numbers — many terms 
of the operator (34) are needed, so that it would be desir- 
able to find other methods of evaluating \]/* for these 
values . 

Call this domain in which the operator (S4) may be ex- 
peditiously employed Ej^. In the remaining part of the in- 
tersection of L and H, xj/ * may be obtained by analytic 
continuation, as well as in that part of S which lies out- 
side L . 

In section 17 of reference 3, the general ideas under- 
lying such a procedure were developed; in the following, a 
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more detailed, discussion of this procedure will he ^iven, 
and some methods indicated for oTercoming certain difficul- 
ties which arise in its application. 


In the following, the particular solutions which are 
due to Chanlygin will he gi-ven, in which case it is conven- 
ient to replace the variable q hy 


(k - 1) ^ 

T = ^ — 2 ^ ^ 

2a^® 

1 

^ 

k - 1 

Obviously, for k = 1.4 and = 1 


T “ 


1 

5 


q®. 


P 


5 

2 


(6S) 
( 69 ) 


( 70) 


the values of T for corresponding M, 3, T, r/a^, \ are 
given in tables 3, and 4. 

If equatiOi. ( 13 ) (from which (I 8 ) was derived) is ex- 
pressed in terms of T and. 6, it assumes the form 




1)1 


2 T (1 ~ t) 




69 


2t( 1 - t) 




5 \|/ 


= 0 


(71) 


Following Chaplygin and separating variables, the solution 
of ( 71 ) may be written in the form 

00 

\|^(T,e) a^j(T)cos + b-p(T ) sin ( 72 ) 

( -L < 9 < L ) 

where a^ , b^ satisfy the differential equation 
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d 

|^T( 1- T ) ^ 

dajj^ 

1 — ( 2p + 1 )t 

(1- 

T)-' 

2 S 
V TT 

dT 

dT J 

T ( 1 - T ) 

41® 

d 

["t ( 1- T )“P 

dhi,3 

1- (26+ 1)t 

(1- 


V®TT®hu 

dT 

dT J 

T( 1 - T) 

4L® 


(>? 3 ) 


which can. easily he transformed into the hypergeometric 
equat ion . 

Unf ortunat ely , the power series for these functions 
converge so slowly that, in most cases, this approach must 
he abandoned for computational purposes. (See sec. 17 of 
.reference 3.) On the other hand, if certain changes in this 
approach are made, numerical results may he obtained. Two 
such procedures will he discussed in the following: 

1. The solutions of (73) are expanded in the form of a 
power series, not, however, around T = 0, which is the 
case of the hypergecmetr ic series, hut around some conven- 
iently chosen value of T which lies inside the interval 
under consideration. ffor example, it will often he conven- 
ient to- introduce, for T, the variable 

CT = T - 0.15 (74) 

Since the particular solutions are often considered only for 
a small range of variation of T, say, 0,13 < 7 < 0,20, 
the coefficients of (73) may he approximated by polynomials 
in a, so that if, say, L = rr/2, the equation approximat- 
ing (73) is for p = 5/2, 

d®U 

(_a® + 0.5500 a 0.2325 a + 0.0191) 

+ ( 1.50000 a® 1.4500 a + 0.1838 ) 

aw 

+ (0.1000 - 6.0000a) v% = 0 (75) 


It follows that every solution of this equation can he rep- 
resented as a linear combination of two independent solu- 
tions, each of which can he written in the' form of an in- 
finite series in a, which will converge for j cf| < 0.15, 
at least, a range of convergence which will he sufficient 
for most purposes. 
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Solutions Ui (ctj v), w) are chosen so that 


Ui(0; v) = 0 
UsCO; u) = 1 


Ui(0; v) = 1 

U^( 0; u) = 0 


( 76 ) 


Thus^ 


Ui(o-i u) = a + 


1 + baO** + b_a® 


+ . . . J 


(77) 


^2(0; v) — J- 'l Ug 

aj, bj^, i = 2, 3, . . . being functions of u . 

If the foregoing expressions are substituted into (75), 
the following sets of eq^uationa are obtained 


0.03826as + 0.18375 =0 

1.45 ~ 0.1 V® + 0.8325as + 0.11-175^3 » O 
1.5+ 6 u®+ (4 - 0.1 u®)aa + 1.94625a3 + 0.?295a4 = 0 


> ( 78) 


0.0l9125(n+l)(n+2)an+jg+ ( n+-l )( 0 . 18375 + 0.235n)an+i 

+ ( 0 . 9 n + 0 . 55n^— 0 .1 u® ) a^ +(— n® +4 . 5n— 3 . 5+6 v® ) = 0 

n = 2, 3. 5, . . ‘ 




^The asymptotic behavior of solutions' and Ug for 

large -values of n ca.n be obtained ea follows: If the coeffi- 

cient y in (18) is a-D-oroximat ed by cX“~ (for k = 1.^, 
c = *035) subsequent formulas become simpler (e.g., see 

Ref. 3t P-, 31 )* particular, equation (Hj) of R“f . 3 becomes 

H®\l/00 + = 0, i.e., 1(E) = (l-M®)/p® = E° , where s is a 

certain constant which is not eoual to 2. For X = 0, i.e., 

M * 1, H = 0. Therefore, if (H) cos nS is 

substituted into this equation, there results 
l(°-) = E:(n®/(®“2^ H) where the functional form of K is 
independent of n. The expression cX~^ ap-oroxlmat es 
asymptotically the exact value of ?, since in the- neighbor- 
hood of X = 0, 


F 


= z 


u=-o 


^-6+2u 
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0.038251(2 - O.lv* = 0 
+ 0.8325b2 + 0.114751(3 = 0 
0.22951)4 + 1.946261(3 + 1 ) 2(4 - O.lv®) = 0 

0.019125C n+l)(n+2)bii+2+(n+l)( 0.18375+0. 235 ?5n)l(jj+i 
+ ( 0.9 + 0. 55n®-0.n;®)l)j^+(-ii®+4.5n-3.5+6v®)bj^„j *= 0 

n. = 1, 2, 3* 5f • • • 




> ( 79) 


Aqu-tions (78) and (79) then give and ae poly- 

nomials in 15®, Since it is often sufficient to consider 
the series (77) in a small interval, say, |a| <0,05, it 
will suffice to employ only a few terms of this series, and 
hence it is necessary to compute only a small number of 
bi. i’ 

2. To employ the second method mentioned, equation 
(73) is written in the form (setting L = tt/2 as before) 


(0.15 + a)(0.85 


. <3. u . du 

a) + (1.5a + 1.225) — 

da® da 


—vS 


/ 0 . 1 - 6a V 

( = 0 

> a + 0 . 15 > 


( 80) 


and two independent solutions Ui(a; u ) , U 2 (a; u ) are deter- 
mined, not by the power series ('?7) but by an approximation 
method. That is, write 


du 

— = V or Au = vAa 
da 


(8I) 


Then 


X 1.5a+ 1.225 > 

0.l-6a N -| 

. V( 0.15+a)( 0.85-a) ^ 

1 V+15 I j U- ' 

' ^ ( 0.15+a)®( 0.86-a) '' J 
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Assume that Act is in a small interval, say 0.001, and 
that for CT = 0 , u = 0 , V = 1 . If (8l) and (82) are em- 
ployed, Au( CT ) and Av(cr) may be determined at ct = 0, 
and hence 


u( O.OOl) = u(0) + Au( 0) 
v(0,00l) = v(0) + Av(0) 


(83) 


from these values of u and v, Au( O.OOl) and Av( O.OOl), 
and u( 0.002), v(0.002) may be found and the entire pro- 
cedure continued until an approximate curve is found for 
u.(cr), 0 <CT < 0.15. In this manner the desired integration 

is performed. 


Both methods described above are appropriate for the 
purpose- of computation, but they are insufficient to give an 
insight into the behavior of the particular solutions when 
U — ^ 00 . 


As has been mentioned, there are instances in which the 
whole flow is subsonic and therefore it is necessary tocon— 
s ider particular solutions only for this range. In this 
case it is useful to consider certain other expressions 
which will be derived in the following, This procedure is 
likewise based on the method of separation of variables. 
However, instead of employing the variables (q,9) and even- 
tually (t, 9) as in the case of equation (7l), the variables 
(^,9) are employed so that it becomes necessary to consider 
equation (l8) once again, Function F has the expansion 
( 26) . 

Let 

vl/* = Uy cos v9 (or Uy sin u9 ) (84) 

where u now has the meaning 


V = ~ , V* = 0,1 ,2, . . . 

(Drop the subscript u to let U = Uy ; U will vary, of 
course, with the choice of v, even though this is not in- 
dicated by the notation.) Thus, obtain for U the ordinary 
differential equation 

U" _ + 4FU = 0 ( 85 ) 
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If S' is written in tiie form given by equation (26), then 
as before, two independent solutions will be determined. 
These will be denoted by Ij( ^ , TjC ® ) ; two cases must be 
distinguished, that is, whether u is or is not an int eger . 
In the latter case. 


where 

,(i) 


= 1 


U' 


(1) V (1) (u+ 3 n)^ 




n=o 


U 


(2) 


( 2 ) ( —u+sn 


dn e 


n=o 




= 1 


\ ( 86 ) 


(1) 


l+u 


( 1 ) _ — Cg— Cidj. 


dp" 


(") 


( ®) 


(s ) = 


1-v 


— Cg— C^di 


( 8 ) 


4+2u 


4— 2u 


( 1) -On~*^n-i'^i “• • •-Cidjj_j^ 


<^n = 


( s) — Cn,— Cn_idx- 


.-C 


i*in— 1 


n +nw 


n®— nu 


Cjj are defined in equation (26) and listed in table 2; 
while in the former: 




n=o 


(3 ) (u+Sn)^ 
e 


U 


(a ) 


•» = rt 




y ( 87 ) 
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The convfrf'ence of. .these series will he di^oussed in annen- 
dix IT. 


In section 17 of reference 3f s- "procedure has heen in- 
dicated for constructing a mioced flow (that is, a flow 
which is partially supersonic) around an obstacle the 
boundary of which is a closed curve. 

If the boundary curve is nrescribed, then, probably^, 
in many instances,, no solution of the problem exists, and, 
therefore, the problem aris'es of finding necessary and suf- 
ficient conditions in order that a (mixed) flow pattern 
around an obstacle exist. In considering this equation, 
two cases have to be distinguished: Zither 

1. The behavior of the flow at infinity is completely 
prescribed (for exam-nle, a uniform flow with no vortlcity), 
or 


2. The flow at infinity possesses certain properties, 
but any solution such that its behavior at infinity is of 
a certain type will be considered as admissible. 

The second problem is to determine when the solution is 
unique^ and, further, under what conditions it is "stable." 
This situation suggests considering a problem Which mathemati- 
cally is much simpler (but still exceedingly difficult), 
and, in certain Instances, can give the type of ,°nswer de- 
sired by the research’ engineer, namely, the quest-ion whether 


^Unfortunately, no definite results, in this direction 
ere known to the author, ' On the other hand, investigations 
of problems of similar type seem to indicate that the above 
is the case. See, for instance, reference (9)- 

■ ®It would be hecessary to- take into account the fact 
that solutions with a free boundary can’ exist. See exam- 
ples of solutions with a free boundary in incompressible 
fluid case given in reference 10. What occurs when no 
continuous solution or a solution with a given free bound- 
ary exists, and under what condition a shock wave arises 
are not known. (See .fig, 6.) 



obstacle. A free boundary is indicated by , 
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to a given hodograph there exists a flow possessing the 
prescribed behavior at infinity. Often it will suffice to 
know not whether a flow pattern exists around some pre- 
scribed (in the physical plane) profile, but only whether 
such a flow exists with' prescribed (or approximately pre- 
scribed) velocity distribution. Then the whole investiga- 
tion may be shifted to the hodograph plane, and the prob- 
lem considered as formulated above. As an example, a flow 
around a Joukowskl profile will be considered in order to 
determine whether to a hodogre-ph similar to that indicated 
in figure 7 there corresponds some mixed flow. In the 
following , s ome necessary conditions for the existence of a 
flow pattern satisfying certain conditions will be given. 



While, in the case of a subsonic flow, it is more 
convenient to operate in the (\,9)— plane (in the pseudo- 
logarithmic plane), in the case of a mixed flow, it is 
convenient to return to the ( H ,0 )— plarie^ , since A(M) be- 
comes imaginary for M > 1. 


As X and H are connected by the relation 

j W 1 ^li S 

dX <= , 1(H) = , it does not present any theoret— 

v^l(H) p® 

ical difficulty to replace the variable X by the variable 
H, in functions S jj.( X , 0 , X q , 6 q ) , k = 1,2,3, which represent 

a pseudo vortex, and two pseudo doublets, respectively. 

Thus, functions are obtained which shall be denoted by 

Tk(H,0 :Ho,0q) = Sk(X,0 :Xo0o) (68) 


^Here. H 


is employed in the sense of equation (32). 


1^0 
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It is now assumed that at the point (Ho,&o) (which is the 
image of the point at infinity in the physical plane) the 
function hehaves like 

AiIiCH.e ;Ho,eo)+Asa?2(H,9 ;Ho,8o)+A3l3(H .0 iHo9o)+R(H,0 ) (89) 

where E(E,9) is a regular function of 

-^s > -^3 must satisfy the conditions indicated in sec- 
tion 14 of reference 3 in order that 3 11 the- physical plane a 
flow around a closed curve is obtained.) How, the question 
arises whether or not a function R( H , 0 ) can be found 
such that the expression (89) vanishes on the boundary line 
of the hodograph. (See reference 5.) 

Remark ; It is not Immedlacely apparent that 9 ; Hq , 9^ ) 

is a single-valued function; however, it is possible to prove 
that this is always the ca.se. The assumption is now made 
that the stream function can be approximated in by 

Chaplygin's solutions 


M 



Ul 
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A system of functions will now be constructed which, as 
will he shown, possesses the property that every function 
which is regular in Si can he appz* oximat ed hy a conven- 
iently chosen combination of these functions. As has been 
proved in appendix I, every solution of (I 8 ) can be repre- 
sented in S 3 _ in the form 

E(H,e ) = Ee[P( g)] 

P( g) = °^g( Z) + Z) + . . Z = \(H)+ie 

(see equation (30)) where g(Z) is a conveniently chosen 
analytic function of a complex variable. Let £1 be the 
image of in the (X,9)— plane, and let fl'(A-l-i0) be that 

function which maps Si into the unit circle (see fig, 8 ) 

W(H,9 ) = W[X(H) + 10] (92) 

maps Si into the unit circle. Every regular ' analyt ic 
function g(X + i9) of a complex variable X. + 10 can be 
represented in Si in the form 

00 

g “ V a 

Li V 
p=o 

Setting 

Esn(H,9) = Ee|p [w'"(\(H) + i9 )]j 
Ji 2 n+i(H, 9 ) = Imj^P [w“(K(H)+i0)]| 

a system of functions is obtained such that every function 
\J/ can be developed (and therefore approximated) in , by 
a series of such functions. 

Eemark ; The classical theorem of Eunge on appr oximat i cn of 
analytic functions states that an analytic function can be 
approximated by a polynomial in a simply connected doma.in, 
while, here, an auproxiraat ion on the boundary also is needed. 

This difficulty can, however, be overcome in a manner 
similar to that described in reference 1?, Since \(/ * setis- 
f-ies in the linear partial differential eaua.tion 



4 ? 
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= 0, the function xI/^CW) in the W-nlane also 
satisfies an equation of the same t j"p e . How, in reference 12, 
it is proved that to every solution which is 

continuous in the closed domain, another solution of. 

the same equation can he determined which is re^rular in the 
closed domain, and such that | - '\1/* j G , where e > 0 

can he chosen arbitrarily small. 

The fact that the domain Hi goes to infinity does 
not cause any difficulty, since the function \|/* for . —X. 
sufficiently large is arbitrarily small and the domain Hi 
can he replaced by a hounded domain, assuming that the 
line — — — — is part of the boundary curve. (See fig. 9.) 



Assume, now, that the reauired function exists. By 
the assumptions which have been made and the foregoing 
considerations E(E;0) may be approximated In Hi by 


N 

Z_* 

U = 1 


a.^^Hy(H,9) and the entire term (39) in Hg 



by Cy ( H , S ) • 


Ther ef or e , 


lim [ 

N -> CO / 

M — > 03 


f . 

■ N 3 ■ 

s 

* 

/ 

y^ a\;Eu(H,0 ) + AyTy 

ds + y 

\ bvCy(H,9) 


_u = i u=x ~ 

■'^hs 

_u=i 


ds 


r N 3 M 

r y auEy(H,e) + ^ ^ byCy(H,9) 

'-'k Lv=i u=i- p=i 


ds 




r U 


U 


y ayRy(H,9) + y AyTy - y byCy(H,0) =0 (94) 


Ln=i 


v=i 


U = 1 



^3 
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Rem ark ! Note that ay = and hy - by^^^ will, in 

general, depend on N, although this is not indicated in 
the notation. 

Rrom the preceding it is then possible to conclude 

t hat : 


If a?^(H,0; Eo,9p), u = 1,2,3 denote functions with 
singularities at (H^,9 q) (see (88)), and Ru(H, 0 ), u = 

1,2,3 particular solutd.ons of (18) which are defined 

by (91) and (93) in the domain L, and, finally, Cy(E,9), 

U = 1,2,3,. • * Ohanlygin solutions which are obtained 

by replacing T in (71) by t = t(B) and annlying the 
method of separation of variables, then, a necessary con- 
dition for the exist"ence of a flow the hodogranh H of 
which is prescribed and which satisfies the hypotheses 
previously indicated is that for M and N sufficiently 
large, and for conveniently chosen ay (= ay^^O and 
h^(= by(^)), the expression on the left-hand side of (9^) 
(omitting the limit signs, of course) can be made arbi- 
trarily small. 


C0N0LUDIN& REMARKS 


The treatment of two-dimensional irrotational motion of 
a compressible fluid, developed in this and preceding publi- 
cations (references 1 through 5), can be considered as a 
direct generalization of the classical methods employed in 
the Incompressible case. 

In the present paper, a representation for a stream 
function of a compressible fluid-flow in terms of an arbi- 
trary analytic function (which had been previously derived) 
has been considerably improved. An analogous formula for 
the stream function of a supersonic flow in terms of two 
arbitrary, twice differentiable- functions of a real variable 
is also obtained. A method is developed for extending a 
subsonic flow defined in a portion of the plane into a 
larger domain. In some instances this process will lead to 
partially supersonic flows. 

The procedures described for determining f 1 oy ' pat t erns 
of a compressible fluid require, as a rule, long coriputat ions 
which necessitate the use of modern computational devices. 


44 
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On th.e other handj the eauation, (24), for the stream 
function is linear, and therefore the principle of super- 
position of solutions may he applied. This fact suggests 
preparing an "atlas" of flow patterns, which should include 
stream functions \J/y(v,0),u = l, 2, . , . ,n of a num- 

ber of basic flows (say, a flow around an obstacle of ellip- 
tical shape, etc.) as well as a number of simple solutions^ 
\{/|j_(v, 6),= u = 1, 2, 3, , . of the compressibility 

m 

equation ( 18 ) . Every combination \j/y + where 

|x— 1 

are constants, represents the stream function of a pos- 
sible flow. By -using some auxiliary tables (see reference 
5) and conveniently changing the ajj^'s, an engineer will 

be able, comparatively quickly, to vary every "basic" flow 
mentioned above, so as to obtain flows approximately, at 
least of the form desired. This procedure can be appl'led 
for supersonic, as well as for mixed, flows. 

Two-dimensional irrotational motion is only a rough 
approximation to the actual situation, since most flows are 
three-dimensional, and friction, turbulence, and so forth, 
influence the motion. 

By the hydraulic hypothesis (see reference 15, p, 84) 
turbulence can, however, in many instances be disregarded. 
Eurther two-dimensional solutions can, in the axially sym— 
meti ic case, be considered as a first approximation and 
used in order to obtain better approximate solutions of the 
three-dimensional problem. By using the two-dimensional 
solution, determining the density p = pQ, and replacing 
p by Pq in the equations VCp"?) =0, V X = 0 (for 
the three-dimensional case) the above system is reduced to 
a linear one, V ( Pq"^) =0, V X q = 0. A solution 

of the linearis^ system in which satisfies the re— . 

quired boundary conditions can be considered as a second 
approximation for the three-dimensional case. 


^In certain cases it would be advisable to use for 
the functions introduced in reference 5. 
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The use of certain tables which need he prepared only 
once will greatly facilitate the work. 

The nonaxially symmetric case may he treated hy com- 
bining the method described aboTe with certain operational 
processes developed by the author for generating three- 
dimensional vectors satisfying 

V "q = 0 and 7 X = 0 

3y employing a procedure similar to the one given 
above, the boundary layer may be taken into account. 


Brown University, 

Providence, R. I., October 1945. 
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APPENDIX I 


PP.OOE OP A PUNDAMINTAL THEOREM ON SUBSONIC PLOWS 


In the previous reports of the author the hasic 
theorem on which the entire method was developed was'i 


Theorem I ; Let P*(2A) he an analytic function of a 
real variable A, defined for ~cs< a ^ A < — c<^» which 
possesses the property that 


d%* 


dA 


Z 


c(K + 1)1 
(e - A )^-"2 


for a<A<— c, K=0, 1, 2, 


. ( 95 ) 


where c is a suitably chosen constant. 

Further, let Q^^^(2A), n = 1, 2, . . . denote a set 

of functions which are defined by the recurrence relations 


( 2n+ + QAA^^^ +4P*Q,^“^ = 0, n = l, 2 . . . (96) 

( 1 ) 


Qa 


~ 4 P*» 


= 0. 


— 00 < a < 0 


Finally, let g( t an analytic function regular in 
a domain B, which contains the origin. Then 


CO 

\l/*CA,e) = Imi'g(z) + y ’ ■ 2A )g ( Z (97) 

L Z_. 2 ''“•n I J 


n=i 




g^°^(z) = g(Z) 


^1 
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will be a solution of 


L(\l/*) = 1 A\l/* + P* \|/* = \i/ * + P* \1/* = 0 (98) 

M* Z Z 

which solution will be defined in every simnly connected 
domain lying in the intersection of B and H, where H 
denotes the domain [ft's X ^ , 6®< 3 X®]. 

This theorem was not, however, "Droved-: for the P of 
(19)* but was proved for Pm. where Pm was a polynomial 
in e^^ 'which approximated P in the interval 
Xq < 0 to a previously specified degree of closeness; and 
it was further shown that these P^j^ satisfied relation 

(95) as well as the other hypotheses of the theorem, 

Tn the follo"Wing, it will be shown that the P of 
(19) actually satisfies the hypothese of the theorem in 
the case where a > -00. 

In order to .prove theorem I for the P of (19)i it Is 
necessary to show that 

Theorem; To every interval 

I - [a = X'^ - €»]' “°°< a,e<0 (99) 

there exists a constant c = c(a, -e)< <» such that the 
function P of (I9) satisfies the inequality (95)* 

The main idea of the proof described below consists in 
showing that P, when continued in the complex domain, 
that is, P(z), Z = Xi ■+• iAg> is an analytic function of 

the complex variable Z, which function is regular in every 
circle 


z + X 


^ - X, 


e 


4g 
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Function F(T) is a rational function of I the only- 
singularity of whioh (for finite values of T) is a pole 
at 0? = 0. Sherefore, in order to prove that F is a regu- 
lar function of Z in the domain^ 


D = E[— o3 < a % Re Z — n < 0, — a>< 2rrk Im Z ^ Suk < oo] 

it suffices to prove that 2 is a regular function. of Z 
in D and does not vanish there. Let X = e". 2hen as Z 
varies in L, X varies in the domain | 

0^e =jxl^r< 1 

with £ = exp(a) and r = exp(— n). 


Instead of considering 

X = ^ ~ ^ / I + h 

1 + 2\1 — h 2/ 



k > 1 ( lOO) 


it will he found convenient to perform the substitution 


s = 1 - 2 


( 101 ) 


All the above-mentioned considerations remain unaltered 
by the change of argument 2 of (lOO) except that it now 
become's necessary to prove that s is distinct from 1 in D. 


By substituting 2 = 1 — s in (lOO) the following for- 
mula for X is obtained: 


.-1 


X = 


'-h-" 


s L h 


+ 1 


— 1 


1 + s J 


( 102 ) 


Unless h ^ is a rational number, X will be an in- 
finitely many valued function of s. In the following dis- 
cussion, only that branch of the function X will be con- 
sidered which has the property that 


Im ■< log 

^ ^ - °i 


0 ( 103 ) 


L h — 1 + s J 

1 

(0 

II 

O 


C 3 denotes the domain defined such that the co- 

ordinates of the points belonging to this domain satisfy the 
conditions given in the interior of the brackets. 
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This "branch is uniquelj' determin-ed within a circle of 
radius h“^ - 1, for the branch point nearest the orifrin 
occurs at s = 1 - h“ ^ . In the following, by X(s) will 
always be denoted the branch of X(s) described above, un- 
less the contrary be specified. 


Since X(s) is regular at the origin it can be devel- 
oped in a power series 


X(s) 


CO 


I 

nal 


anS 


n 


ai 0 


( 104 ) 


H"ow, as ai ^ 0, the inverse function 
X = 0 and can be represented there in 
series 

CO 

9 (X) = b^X^ 

n=i 


s(X) is regular at 
the form of a nower 


(105) 


It will be shown that sCx) is a regular function of 
X for |x| < 1 and therefure by a classical theorem of 
analysis, the power series (I05) converges for |X|< 1 , 


Lemma ; The function s = s(X), which is the inverse 
of the branch of the function X = X(s) of ( 102 ), which 
satisfies condition (103). is regular in the unit circle. 

Therefore, at" every point inside the circle Isj < 1 
the power series can be inverted and js | expanded as a 
power series in X - X^ , where X^ = X(sq): 


s = Sg + bi(s^)(X - Xp) + . . . (106) 

It is now necessary to prove that the Inversion about 
8 = 0 converges throughout the int’erior of the circle 
I X[ = 1 , Two independent proofs of this fact will be pre- 
s ent ed. 
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A. Determine the values of q at which the inversion 

of X(a) may hecome imposaihle. These points are at s = .1 1 

2, -1/h, + 1, i + 1, and s = co, the corresponding 

h h 

values of |X| hoing 1, oo,os, 0, and 1, respectively. 

Since only that branch of s(X) for which s(0) = 0, not 

s(0) = 1 + i , is to be considered, it follows that the 
h 

function s(X) determined by the inversion of the specified 
branch X(s) is regular inside the circle |X[ = 1. There- 
fore the expansion of s as a power series in X will 
converge for jxj < 1, for as has just been shown by the 
foregoing reasoning, for no value of X Inside the unit 
circle can s(X) be singular. 

B. Clearly, X(s) = 0 at s = 0 but nowhere else in- 
side or on the circle |s|- 1. In particular, on the cir- 
cle I s| =1 it is not difficult to chow that the minimum 
value of |X| is 1. Hence, by the theorem on page 136 of 
reference 13 and the proof presented there it follows that, 
since dX/ds does not vanish at s = 0, s can be expanded 
as a power series in X convergent within the circle 

|Xt =1. 

z 

Since X = e and since it has been shown that s is 
an analytic function of X in |X| < 1, it follows that 
8 is an analytic function of 2 for Re(z) < 0. 

Finally, to prove that s cannot be eq^ual to 1 in D 
(and hencer that T 5^ o) , assume thf contrary: namely, 

s = 1. Then by equation (l), |Xl = 1, and therefore 
Re(z) = 0. But, for all Z in D the real part of Z is 
negative. Therefore T is distinct from zero throughout 
D, and thus F is a regular function of Z in D. 

As has been indicated in section II, the proof of the 
validity of the representation (24) for |0j < X < 0, 

has been given under the assumption that the quantity a, 
for which Q^’^^(a) = 0, is larger than -co. 

In the following the proof for the excluded case a = —00 
presented. The proof is essentially based on the 
following corollary to the theorem, equation (99). 

Oor ollary ; To every positive number p < 1 and every 

negative numbei there exists a number c(p,x(o)) such 

t hat 
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K, 

d?\' 


d 7 
K 


' (-y^v) ^ ® 


K = 0,1,2, . 


Pr 0 of ; Let G(z) = / P(^)d^,Z = ^J^+i\g. As has 


"been shown above (see also pt . II of sec. 15 of reference 3 
and appendix I of reference 5), P(Z) can he represented for 
Xj < 0 in the form of the following series! 


PCZ) = 


CO 



a. 


„ snZ 
e 


n= 1 


(The constant term vanishes since F(-oo) = O.) Therefore 
(justifying term-hy-term integration hy the usual argument 
involving uniform convergence) G(z) may he expressed in the 
form! . . 

os 

S(Z) = X 

nai 


Having chosen it follows from the last equation that 

' ( 0 ) 

there exists a constant, A(X ) such that for 


Xi = EeZ < X 


(o) 


0(Z) 


< A(X^°^)e®^^ 


Now consider any real X such that X < X^°^/(l >- p). 

(It is understood that X^°^ and p are held fixed during 
the present discussion.) Draw a circle c with center at 
Z ■= X and with radius (-pX); this circle lies entirely to 
the left of the line X^^ = X ° . According to the Cauchy 
integral formula; 


6 la 
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<iZ^/z=A 


iZ I Z=x 


(g + 1)1 
3Tri 


G(Z)dZ 


(Z - A) 


K + 3 


By setting Z = A - pAe^^, 0 = $ = 2rr, the last equation 
hecome s : 


an 


d%\ ^ (g + 1) I p (.-pA)G(A - pXe^^)e^^ 


^ ,rg I 

,<^Z /z = A 


2tt 


J 


(-PAS^^) 


K + s 


N ovr : 




= -PX, 


and 

(-pA)G(A - pA e^^)| < 
Since, for 




0 < <f> < 2tt. He(A - pAe^®) < A(l - p), 
there is obtained the inequalitys 

|(_pA)G(A - pAe^®)| < (-pA)A(A^°^e®^^l-P^ 

Since A < A^°^/(l - p), it is apparent that there exists a 

constant c(p,A^°^) such that the right-hand side, of the last 
inequality is less than c for all A. Therefore, from the 
second form of the Cauchy formula given above, there results; 


K 

d P 


dZ 


K 


Z=A 


< c(p.X^°^)(g + l)i 
(-PA)^'"^ 


A < 


1 - p 


which yields the statement of the corollary. 

Now, given any A < 0, let A^°^ be chosen as follows: 

A^o) <0 and p positive but less than one such that 

A < A^°^/(l - p), and then the corresponding number c(p,X^°^ 
has to be determined. 
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Further, the dominants have to he defined hy 


= 4c /' - 


4c 


(~PA)^ P^(“X) 


and the recursion formula 


• • • 


A 

(2n+l)Q^'^J;'’^ = + 4c(-Xp)‘^ y dX, «) = 0, n « 1,2, 


rather than hy"^ (94) of reference 3. 

Now it will he shown that there exist a set of constants 


.(n) 


such that: 


~(n) (n)'/ .-(n+i) 

Qx = c (-pX) 


(n * li2f« • .) 


Proof : This equation is seen to hold for n » 1 with 
0^^^ = 4c. Suppose the theorem holds for n = jiq . Then: 

X 


— 00 


By carrying out the integration, there results: 

(2no + = (_pX)’^o+e ^ + i) + ° 

By letting 


(no) 


P II, 


t ^ (hq)/ (hq)/ 

(uq+i) p c (Uq + 1) + 4c c /p Uq 

(2no + 1) 

it is seen that the induction is complete, and the statement 
made above is estahlishod. 


^It should he remarked that in formula (94) there is a 
misprint : (—X) ^ should he replaced hy (— a)“^« 
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If the last equation on hoth sides is divided hy 
the following fornula is obtained for the ratio of successive 
coef fici ent s : 


" 7 ^ 


p (np + 1) + 4c/p np 

Sno + 1 


By letting ng — ^ «, the result obtained is (dropping the 
subscript on 


Lim 
n — oo 


(n+» ) 



V n~(n ) , . 

1 + t ^ 

. i ■ n = 1 


from which the convergence of the aeries 
(see ( 97 ) of reference 3) is assured for 
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APPEND IX 1 1 

THE EVALUATION OE THE Q,^ ^ ^ ' a AND THEIR DERIVATIVES 

In appendix I, it has “been shown that operator (24-), 
when applied to a function g which is analytic in a 
domain B, generates a solution of (iS), which is defined 
in that simply connected subdomain of B, the coordinates 
of which satisfy the relation 


0 < A < 0 

In this fashion, possible stream functions in the (A,0)— 
plane, of a compressible fluid have been formed, 

Hovrevor , it is necessary to evaluate the functions 
^ ^ , defined by (lOV), in order to apply operator (24) 
in pract-ice, and in reference 3 this has been done for 
Q^( 1 ^ q(*')^ and . Often, a greater number of 

0,' s is required, and for this reason, in this appendix 
q_(n), n = 1, 2, . . 4 , 8 are computed. 

As it; naa jeen shown on pp . -I*’ s-s defined, by 

(ig) does actually satisfy the conditions of theorem I of 
appendix I, when a = -<», The can be determined as 

follows: 

( 2n+i)Q;^^“'^’ ^ + 4E = 0 (107) 



—00 


Q,^°’^(— o=) = 0, n = 1,2,. . , 

If equation (20) be differentiated, it may be easily 
seen that 


d\ ' -5T^ 

dT ( TS_i )( -ps_6 ) 


( 108) 


so that if E be given its value in equation (19), then 
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—4 


p cLX = — 4 


F — cLT 
dO? 


X=— c 


T*i 


-4 j {-0.12T^+0.51+0,P.1T ®-0.6SI~* + 0.455?“^ ) 
T=l .. / -5!r® 


(■ 


) 


dT 


( )( T ^-6 ) 

-0.50T"^ +0.35T“=^-2.40!D+1 .44 + 1.28 logf' ^ (l09)^ 

^^ 6-1 




(i) 


iNote that there is a misprint in formula (107) for 


giTen in reference for general k, the correct 
formulas for and 4'® ' ere: 




Cl) = Cl+k) 

a 


( l-3k) 
I- k-1 


j +. 5> 1 __ ^5 ® ^ 

k+1 3 h(k+l)®(k-l) 


X log 


hT+1 
-hT+1 J 


T 


1=1 


( 109a) 


^C® ) _ _ 3- ( k+1 ) 

2048 (k-l) 


(3k-l)' 


j 3 


( 6k®+28k^»->58k+16) 
(k-l) 


( 41k^+96k® + 18k?-112k-91 ) (k-l) 

(k+l)=" 

„ ( 44k°+193kg+336k+80)( k-+l ) ^_3_ ( 39 k^-80k-115 ) ( k-l) ^_s 


3( k+l) = 


5( k+l) 


PR 

- lOk(k-l) — (k®-l)T“® 

9 


64h 

( k-l)“( k+l)' 


log 


l/h-T 

l/h+E 




( 109h) 


T =1 


5S 
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B'rora ( 107 ), for n = 1, it may be seen that 





+ 4 


dX' 


_4 II _ q^ 



)dX=/ d(4F.i4'*)*) 


dX 


dX 


—00 

Th er ef ore 


q^ ^ ^ = i 45- + iq^ ^ ^ = • 43 , ^ iq( ^ ^ 

3 V 2 / 3 6 

“CO 

In ( 107 ), set n = 2 to obtain 


5q[® ^ 


( 110 ) 


(111) 


( 112 ) 


and if this be integrated frpm --<» to X, the following 
expressions are obtained 

.X 

5 q(s ) = _q^ - 4 


rq^ ®^dX 


= -q 


( 8 .) 16 


jsdx + i r qs ^ 


ij. 


dX 


dX 


( 113 ) 


If q^ ^ ^ is given its value in ( 109 ) , (113) upon in- 
tegrating the last term assumes the form 

T 

(=) = - 




90 


15 


dT 


dT 


T=i 


= -D.2000q[®^ + 0.011l(q^ f +0.0256T®—0.1152T+0 . 219 4 + 0 . 04531—1 


-0,15 750? ®-(-0.03 760?~®+0. 04200?’"’’-0.02000?“®-f0.08 8 6 
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R emarlc ; and Q. ^ ^ ^ can also "be expressed in closed 

form as functions of T . 

If n is set equal to 5 in (l07), then 

.(4) ^ J3) 


= -q: + 


(-4F)<i^^^dX 


and if (—43’) is replaced "by ^ ^ and the indicated inti 
gration by parts performed* then the formula 


70,(4 ) - _q(^3 ) + q(i)q( 3) _ / q(i)Q(3)dA (115) 


?^=—oo 


Cz ) 


is obtained. Now, insert for the expression 


0,(3) cs — — .i pq( ®) 

\ 5 W 5 

obtained by setting n = 3 in (107). This yields 

7q(4 ) ^ _qU ) ^ (j(i )q(3 ) ^ I r q( d\ 

^ 7 \ =—00 


+ ly 


( 116) 


In order to evaluate the first integral, two successive 
partial integrations are pez’formed: namely, 

= + J* (117) 
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Insert this expression into (116) and employ (l07) for 
n = 1 : namely, — to obtain 

7 ^( 4 ) = I ^ + |- rq^®^ 


+ ^ 4rq^^^)dA 

= q^^) q(^)q^=^>4 i q(^>qr^^ ^ ^ (US) 


Eemar k i By using the previous results, (II 8 ) may be ex- 
pressed in a closed form as a function of T. 

It has not as yet proved possible to express Q/ \ 
n > 5, only in terms of the preceding q(^'*s and their 

derivatives, which formulas would yield for q( ^ ® » 
closed expressions in 2, Cone eouent ly , it has been neces- 
sary to derive formulas for q(^^» n = 6 , 6 , 7, 8 in terms 

of the preceding q^^^’a and their derivatives but which 
also include one integration, 2 hus , a formal computation 
leads to the following 


T 


<q 


( 5 ) ^ _ 1 0 ( 4 ) _ 4 
9 ^X 9 




( 119 ) 


q 


rn 

(®^ = — ^q^^^ +— q^^^q^®^ — ^ /*q( ’■^q(®^— dC (120) 
11 ^ 11 11 J dT 


q( 7 ) ^ -q( + q( ^^q( — q( - — q^ q^ 

^ 11 ^ 22 


+ 




(^^5 ) q(s ) 


d2 


dT 


( 121) 
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15 


195 


“ ~— ^ ^ ^ Q.-.^ ^ 0 ^ ® ^ 


390 


Q. 
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dT 


dT 


( 122 ) 


In. order to evaluate the expressions obtained for the 

Q^^^’s, not only the pracodine ^ Must he known as 

functions of T, hut their derivatives and the derivatives 
of S' must he known as well, 

The following tables supply some of the needed deriv- 
atives. Observe that 


0 . 


for S = 0, 


( m ) 


m = 1,2, . 

S >= 0,1, . 


( 123) 


0 , 


In table 5, typical functions necessary for the ex- 
plicit evaluation of the Q's have been tabulated. 
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Table 

Derivatives of of order 


^(1) 


»(S) 


q(3) 


S 

0 

s 


-4^ FdX 
^(s-l)p 


-4 


2K 


(s-l) 


F + 


( 1 ) 


+ I F^- I FQ 

4 F - 4 F qC^) + 16 F® 

3 n2 3 -S' 

4 F , - i F + 16FsF 

3 \3 3 A 

I - I t M + 16(F^)^ 

4 ^ 4 -p. n(l) + ^ F F + '^v- F r>F 

t '■^6 - t ’■>>5'*^^’* * W^\S>‘ 

- i ■" A'**'”’® - 


1 ,(|) - I ,(®)f 


- i Q<"’ 

■5 


4 ,r^F - I 


(S) 

- ^ 
5 
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,( 5 ) 


)(4) 


,(5) 


,( 6 ) 


>(7) 


,(S) 


1 i 

✓ A 


- 1 q^ 3)+ a<2^ . ^ ^ 


1 <j(|>- I fq(2) 


7 - 7 

5 


4 fl 


s-1 




I 

- 1 ^ (Q^^^F) 

■ II n J' 

— 0^ 

1 — Q^®^F 

II ^A2 11 ^ 

- 1 q(5) - 4 3°'^ (q(5)f> 

n ^;,s+i n 

tf\ 

-A «'=>(«(»)’*. its 

’CO 

1 - -i 

IS ^X2 13 ^ 


F 

s-1 


1 n(®) _ 4 _3 (Q^®^F) 

^ A 

^ ,W* i ,w,w. ^ ,Wq«- 


_ 1 q(7) _ 4 

15 \s+l IS 
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Table 

The Values of F as a Function of T 


-0.12T® + 0.51 - 

j2 j6 

0.048T® - 0.336T + 

T m3 mb t” T° 


-0.0288T^ + 0.2688T^ - 0.6024 + P » 3.680 + g,., 4,560 

^2 qi4 


- S4.171 93.895 _ 97.978 ^ 34.992 


^6 


t8 


jlO 




+0.0230T® - 0.2688T® + 0.8909T - P-- + S . 2944 __ 59..9,5fi 

T t3 t5 

+ 453.86 ^ 1493,6 . 2357.1 _ 1763.6 . 505.88 


n9 


ill 


nl3 


a5 


■0.0230T® + 0.32261“^ - 1.4454t2 + 2.0993 + 

l2 

_ 70.286 ^ 1063.4 _ 7496.1 ^ 27817 _ 57015 
l4 I® I® T^P 


+ 64722 _ 58094 ^ 9069.9 
ll4 ll6 ll8 


(0.0276t'^ - 0.4515T® + 2.5503T® - 5.5954T + 3.4356T”^ 
-59.495T"® + 1672. 4T“5 - 21263T“'^ + 147246T"® - 598241T“^1 
+1472885T“^® - 2211474T~^® + 1973284“^"^ - 959960T"^P 
+195910T~®^) 


(-0.0587T® + 0.7225T® - 4.9235T^ + 14.540T® _ 0.2690T 
-17.522 + 3.7520T ^ -30.572T ^ - 16.348t"'^ + 1909. 9T~^ 
+22.512T-® - 41666T“® - 9.6480T“'^ + 483449T"® 

-3350057T + 14632715T~^^ - 413377621“^^ + 761271681"^® 

-90418589T~^® + 66612802T-^P - 276468701“^® + 4936939T"24) 
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APPEND IX 1 1 1 

PROOP OP A PUNDAMENTAL THEOREM ON SUPERSONIC PLOWS; 
EVALUATION OP THE E^^^ 


The following theorem was enunciated without proof in 
s e ct i on I II; 


Let 





0 


(48) 


where n(^) (fj is considered as continued for complex 
values of the arguments) is an analytic function of the 
complex variable J , which function is supposed regular 
for [^[ < Ai then, there exists a set of functions 


, n = 1,2, . . .A = (49) 


such that 

P(l ,Ti) = Vid ,n) + Vgd ,n) 
is a solution of (48) where 


(50) 


00 

Vid ,d) = f(i) + y ) 


n=i 


> (51) 


00 

v^U.Ti) = e(n) + yE^UA)g’^’=^ (d) 


and 

f*^°^U) = f(|) 


n= 1 


gl^®3(^) = g(T,) 


f[a+i](|) = / fM(|jdlx. g[n+i] (n) = / gM(ni)dTi, 
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f and g iDeing two arbitrary, twice dif f er ent laTsl e func- 
tions of I and T) , respectively. 


The following is a proof of this theorem. Substitute 
(50) and (5l) and (5?,) into (48) to obtain 

r 

\ j- ■f’fc >ij'. v—yj-rix n \ j I j _ 

&• 


A) 


|^g( p ) + f (I ) 

r 

^(^n+i ) ^ jj(^n^) ^ \ a ) 


^ + E +Q( A ) ^ ^ A ) 


gfi](p) + f[i3U ) 

. . , + (p) + (t) 


B*' 


= 0 


( 124) 


1 (m) 


Thus, if the ra - 1, 2, , . * are determineA by the, 

recurrence relation 

A 


^^(A) = - / a(Ai)dAi 


A 


> ( 125) 


E 


( n+i ) 


,(n) 


(n), 


(A) = - / ^(Ai)aAi 


■'o 


J 


Then (50) v/ill formally satisfy (48). The uniform con- 
vergence of this series will be proved by the method of 
dominants . 


Lemma ; Let Q( ^) be an analytic function of the com- 
plex variable v/hich is regular in the circle l^j < A^ ; 

then, 




< 


2MAq ^ 2MAq^ 

Ao-uriAo- f 


d“n( 0 

"Af' 


^ 2n'.M Ao 


2( n+l) IMAg® 
- , I n+ a 

Aq *" ^ I 


>0 


< S< ^ 0 , i . e . , 

t, real and 
pos i tive 


wh er e 


( 126) 


( 127) 


M = max |f( £)| 

The proof of this lemma follows 
f ormula 


for |t| = A^ (128) 

immediately from Cauchy's 
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n(t ) 

and the fact that 


1 r 

2TTi J i i 
Ul= Ao 


( 129) 


Ao 

Ao-t 


> 1 


for t real 
p o e i t iv e 


and 


( 130) 


In the following O will he considered only for real 
non— negative values of that is, O will he taken as 

a function of Ai 0 < A < Ac* Ihis restriction on A will 
he understood for the remainder of the proof unless the 
contrary he stated. 


Let 


n(A) 


_2H^ 

(Ao-A)^ 


( 131) 


Equat i ons ( 126 ) , 

form 


(iP?) may then he written in the 


|n(A) < 5(A) 


d'^(A) 


d^(A) 


dA 


n 


dA 


n 


( 132) 
( 133) 


respectively. By definition, if (l3p) and (l33) hold, Q 
‘s a dominant of which fact will he symbolized hy 

» Q or n « Q- 


'"( -l) 

If 1^-^HA) is given hy 


^ ^(A) = / n(Ai )dAi 

■u 


(134) 


Ihen 


E^^^(A) 


A 


n(Ai)dAi 


5(Ai)dAi< b‘'^^(A) (135) 


and also 
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cL^e( 

dA^ 




dA' 


Thus 


^^(A)« ^^(A) 


Suppose, now, that 

,A 


( n+i ) 


(A) = 


i|j^(Ai) + Q(A)i^“^(Ai) dAi + hI’^^A) 


where 


jj( n ) <.< g( n ) 


»“( n ) 

0 « ‘ 


Then, it follows immediately that 


n+ 1 ) 


(A) 


< E 


( n+-i ) 


(A) 


dj.( n+i )( 


dA 


dE^ n+i )( 


dA 


and by considering the corresponding derivatives of 

E^^^ +f2E^^^ in comparison with it 

lows that; 


e( « E 


AA 

( n+i ) 


which completes the proof by induction. 
If, now, is given by 


2. . I- 


^in) _ CnAp M 
A 


H 


LAq •” A Ao "* 


» 0 


( 136 ) 

( 137) 

( 138) 

( 139 ) 
( 140) 

(141) 

(14S) 

fol- 
( 143) 

( 144) 
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where is some cottT en isntly chosen positive constant 

(to he det ermined, later ) . Then an explicit expression may 
he obtained for 


,A 

^^CA) = / fi( Ai)dAi = 


•’(A) - 


■Ao - A Ao-1 ° 


(Ao - A)“ An’^J 


> (145) 


It is then seen that 


Cl = Ao® ^ 


(146) 


to obtain the recurrence relations holding between the c^, 
write 

~(n+l ). . N ■ r ( (n+1 )n 


E 


(A) = / ,0 


n 


o 


\(Ao-Ai) 


_^_loUL_ 

(Ao-Ai)^"*® 


dAi + H^“^(A) (147) 


and employ ( 144) , so that 


°n+i ” 


CAo - 1) Ao = 


[ n( n+1 ) + A 1 M] 


( 148) 


( n+1 ) 

Note that 

Cn < M*( n+l) ( 149 ) 

where M* is some conveniently chosen constant. 

To complete the proof for the uniform convergence of 
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'f(|) is assumed to "be a differentiable function 
and therefore there exists a constant, say, Mg, such 
that 



fU) 

< Mg , for 1 

^ 

( 150) 

then 


< for 

— n 1 

CQ 

V 

( 151) 


Consequently, from (5l) it may ha seen that 

/6 „ r, 2 UI M‘^21 ^ 2\i\ ^UHn+l)\ + . . 

Vi(| ,ri)« Mg 1 + + . . . + (15?) 

^ (Ao - A) >i'*(Ao - A)^ -* 


which will converge if 



< 1 


( 153) 


so that the series ’®'i(^,Ti) converges uniformly in the 
domain (153). She same may he said for the series 

hence the series for ?x(5*'n) icay he differen- 
tiated termwise and therefore the formal solution (50) is 
actually a solution of (48) in the domain ( 153 ) . 


Che same also holds for Vg(|,n), hut here the domain 
( 153 ) has to he replaced hy 


ijOj 

|Ao - A{ 


< 1 


( 154) 


Iherefore V(|,'n) represents a solution of (48) in the in- 
tersection of ( 153^ and (l54). (See fig. 10.) As has been 
pointed out in section III, this theorem cannot he applied 
direutly to equation (4?) as given hy (45) has a pole 

for A = 0. In some cases, however, it is possible to over- 
come this difficulty hy shifting the origin. Let a he a 
positive number and 

I* = I - a/? 

T[* = t] ~ a/? 

so that A* will mean 


( 155) 
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S^lgura 10. 


A * * A — a 

Equation ( 47 ) then aasu»e« the form 


(156) 


-i- V » 0 

wh ere 

+ 1, n* + f) 

Fs(A*) = Fi(A* + a) 
so that jPg is analytic for A* = 0« 

The corresponding changes in E^ ^^(A) will he 
cated hy writing ^^^^(A), thus 


eO(A*) = - / B-g(Ai)dAi = - J Ei(AJdAi = Ei(A) 

_g(n+i)(^^.) ^ 


^ TT.Cn+i), 
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Ihe domain of conrergence is merely shifted. 

It is desirable to have explicit formulas for the 

and their derivatives, as these will he employed in 
computation; if n is given the values n = 1,3, 3, 4, in 
( 159 ) , then 

A 

A) = - / 2’i(Ai)dAi (160) 

*“a 

= ?i(A) - I’jCa) + (161) 


e( = )(a) - ri(A) . ilAtl 


6A 


5A 


- f Ei®(Ax)dAi-Ei(a)E^"^A) + i(E^'\A)) = 

' a 

=-(2i!2iA> e^'>(a) e< = >(a) 

N ftA ^A / 


( 162) 


+ E^'-^A) — E X 3( A ) - E 1 ="( a ) ) 

0 1\ V 


2 - - 8 “ 


wher e 


+ ~(E^^\ A))®EiC'A) 

O ' 

9Ei(A) 

dA 


= -Ej(A) 


.E<=>(A) _ 


a A 


9A 


( 163) 
( 164) 
( 165) 
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^(A) 

SA 


Ha) 

aA 


Expl icit 
t i on of 
B( a ) is 

e^^\a) 


aPi(A) 

6A 




SFi(A)l^ ^\a) 


aA 


- 5-i®(A) 


U66) 


= - F,(A)E<=)(A) . . S^^ . .(A)lt . ^> (A) 

aA=^ &A® 


SEi(A)E^^\a) 3aEi(A)’ 


dA 


aA 


ai*^ ^^(A ) 

Ei(A) + Ei (167) 


aA 


formulas for some of these quantities as a func— 
B are given helow; is defined in (45); 

the value of B corresponding to A = a. 


k+1 


16 L3 


'5 b- 3 ^ 2k 


( fc+l) 


B 


(k-l) 


8 


1 B 


h( k+l)=®( k-l) 


tan 


-.1 


hB 


B(a) 


( 168) 


_ f30( fc+l)B“'^+48kB~®+2( 6k-14)B“3 + 2( 3k-l)B 

^ 128 -^ L 


X 


“( l+h®B®)( 1+B® ) ■ 
B=(l-h®) 


( 169) 
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1 ^ ( A X ) d A 1 — 


a 

A 


k+l _ 25 
1024 L 9 


-(k+l)® - lOk(k+l)B 


-7 


( 39k^-80k-15 ) t,-5 . 4( llfc®+4Bk"+59k+20 ) ^,-3 

i ® iTJTIl 

(91-70k-151k2M-52k= + l33k^-14k'^-41k® ) _i 

■ I ■ . S 

( k+l)2( k-l)® 

( l-3k)( 16+6k-12k®-2k®) 

(k-l)« 

(k+l)(l-3k)2 ^3 


B 


3( k-l) 


+ 


128 


h( k+l)®( k-l) 


tan“^ hB 


nB 

B( a ) 


( 170) 


Remark Sometimes the range of variability of the speed is 
comparatively small. In these instances it is useful to , 
replace if' by a constant, say 3’^. Equation (43) then 
becomes 


d®T 

aUri 


Fo? 


0 


and its solutions can be written in 


the form 


t = i 

V(|.Ti) ^ I 

t=---i 


cos 


2t( In ) 


1 /s 


F 


]} H 


F„ 2 


+ e 


rri (i-t®)l 






dt 




where f and g .are two arbitrary, twice continuously 
differentiable functions of one variable. 
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APPENDIX IV 


CONVERGENCE OP IHE SERIES OBTAINED AS SOLUTIONS OP (SS) 


The convergence of the series (86) or (87), which have 
been obtained as solutions of (85) may be demonstrated as 
follows : 


Let 


e®^ = 2 


( 171) 


the half plane Re(X.) < 0 then corresponds to the circle 
I z| < 1 . ■ 


Since 


and 


dU _ dU dz _ ^ eA ££ _ 

dX. d 2 dX dz dz 


d®U dz d sX <3. 

dX*. dX dz ^dX/ 


<i / iiU\ / 

— (2 z ) = 4 z ( ! 

<=l T X X 


\ 


?■ (172) 


d®U dUv 

— *i;) 




( 


dz ^ dz' ^ dz 

the differential eouation (85) may be written in the form 


d^U dU 

+ z — + 

d z ® dz 




U 


( 17S) 


Since P — is a polynomial .in z = e®^ , the point 

z = 0 is a regular singular point of. the differential equa- 
tion ( 173 ) , See reference 14, where it is also shown (in 
sec. 10,13) that the general solution of ( 173 ) may be ex- 
pressed in the following form, provided ' p' is not an in— 
t eger . 


U — CxUi + cb"^s 

where c^ and cg are arbitrary constants and 
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Ui ;= + y (175) 

n=i 

Us = z'^^^fl + y 6nz’^j (176) 


where the ct.^ , are properly chosen constapts. Since 

the only (finite) singular point of eq.uation (l73) is at 
z=0 , the theory presented in reference 14, section 10.13, 
shows that the series converges for all values of z. 

However, if u is an integer, as in the case under 
consideration, the above-mentioned method fails; the ser- 
ies (174) may be retained, but, as shown in section 10,15 
of reference 14, the series (176) must be replaced by 


■’■(I 


\ 2 


n~v 


n-v 




00 

I 


n=o 


n=u+i 


n -u J 


( 177) 


where the are properly chosen constants. Therefore, 

in the case under consideration, the general solution of 
(l73) valid for all z, will be of the f^orm; 


U = CiUi + CsUg“^ (178) 

It is seen that, by replacing z by , (l75) and 

(177) will assume the forms given in (87); the factor X 
in the first term of the second series of (8?) arises from 
the logarithmic term in (l77). 

Since the expansions (l75) and (l7?) are valid for all 
z, in particular for j zj < 1, the expansion (87) will be 

valid for |e®^l < 1; that is, Se(l) < 0, 
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Table 1 


1 2 A 

T and T as functions of e 


oo 


OO 


If T = ^ ^ ” YI. 


2nA 


n=0 


n=n 


then 


n. 

-a^ 

bin 

n 

0 

-1.0000 

1.0000 

1 

0.2392 

0.2392 

2 

0.1087 

0.1659 

, 3 

0.0658 

0.1315 

4 

0.0456 

0.1108 

5 

0-0342 

0.0968 

6 

0.0270 

0.0865 

7 

0.0220 

0.0786 

8 

0.0185 

0.0724 

3 

0.0158 

0.0672 

10 

0.0138 

0.0629 


Table S 


The Coefficients of the Series Expansion 
of F in Powers of e^^ 


F 



Cn® 


2n^ 


Cq = 0.0000 
= 0.0000 
Cg = 0.1375 
C 3 = 0.S858 
C 4 = 0.4335 
Cg = 0.6073 
Cg = 0.7241 
C^ = 0.8678 
Cg = 1.011 
Cg = 1.153 
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Table 3 (cont'd) 


Table 3 

The corresponding values of M, T, "i , q/a^, A 
for subsonic values of II 


M 

T 

T 

9/a^ 

■I 

M 

T 

r 


Iillll9i 

0.00 

1.0000 

0.0000 

0.0000 

+ OO 

0.75 

0.6615 

0.1011 

0.7111 

0.1171 

0.05 

0.996B 

0.0005 

0.0500 

2.6282 

0.76 

0.6499 

0.1036 

0.7196 

0.1093 

0.10 

0.9950 

0.0020 

0.1000 

1.9376 

0,77 

0.6SB0 

0.1060 

0.7280 

0.1018 

.0.16 

0.9887 

0.0045 

0.1497 

1.5364 

0.78 

0.6258 

0.1085 

0,7365 

0.0944 

0.90 

0.9798 

0.0079 

0,1993 

1.2548 

0.79 

0.6131 

0.U20 

0.7448 

0.0875 

0.26 

0.9683 

0.0124 

0.2485 

1.0396 

0.80 

0.6000 

0.1135 

0.7532 

0.0807 

O.SO 

0.9539 

0.0177 

0.2973 

0.8667 

0.81 

0.5864 

0.1160 

0.7616 

0.0742 

0.35 

0.9368 

0.0239 

0.3458 

0.7240 

0.82 

0.5724 

0.1166 

0.7699 

0.0680 

0.40 

0.9165 

0.0311 

0.3938 

0.6036 

0.83 

0.5578 

0.1211 

0.7718 

0.0618 

0.45 

0.8930 

0.0390 

0.4416 

0.5006 

0.84 

0.5426 

0.1237 

0,7864 

0.0560 

0.50 

0.8660 

0.0476 

0.4880 

0.4120 

0.85 

0.5268 

0.1263 

0.7945 

0.0505 

0.51 

0.8602 

0.0496 

0.4972 

0.3956 

0.86 

0.5103 

0.1289 

0.8027 

0.0452 

0.52 

0.8542 

0.0513 

0.5065 

0.3798 

0.87 

0.4931 

0.1S15 

0.8108 

0.0402 

0.53 

0.8480 

0.0538 

0.5168 

0.3644 

0.88 

0.4750 

0.1341 

0.8189 

0.0354 

0.54 

0.8417 

0.0551 

0.5249 

0.3495 

0.89 

0.4560 

0.1368 

0.8269 

0.0328 

0.55 

0.8358 

p.0569 

0.5336 

0.3560 

0.00 

0.4359 

0.1394 

0.8349 

0.0266 

0.56 

0.8285 

0.0592 

0.5432 

0.3208 

0.91 

0.4146 

0.1421 

0.8429 

0,0225 

0.57 

0.8816 

0.0610 

0.5524 

0.3071 

0.92 

0.3919 

0.1448 

0.8507 

0.0188 

0.58 

0.8146 

0.0630 

0.5614 

0.2938 

0.93 

0.3676 

0.1475 

0.8587 

0.0153 

0.59 

0.0074 

0.0651 

0.5705 

0.2808 

0.94 

0.3412 

0.1502 

0.6666 

0.0120 

0.60 

CT.aooo 

0.0672 

0.5795 

0.2682 

0.95 

0.3123 

0.1530 

0.8744 

0.0091 

0.61 

0.7924 

0.0693 

0.5885 

0.2560 

0,96 

0.2600 

0.1556 

0.8821 

0.0065 

0.62 

0.7846 

0.0714 

0.5975 

0.2440 

0.97 

0.2431 

0.1584 

0.8892 

0.0042 

0.63 

0,7766 

0.0735 

0.6064 

0.2325 

0.98 

0.1990 

0.1611 

0.8976 

0.0023 

0.64 

0.7684 

0.0757 

0.6153 

0.2212 

0.99 

0.1411 

0.1639 

0.9053 

0.0006 

0.65 

0.7599 

0.0779 

0.6242 

0.2098 

1.00 

0.0000 

0.1667 

0.9129 

0.0000 






0.66 

0.7513 

0.0801 

0.6330 

0.1996 






0.67 

0.7423 

0.0824 

0.6418 

0.1892 






0.68 

0.7338 

0.0847 

0.6506 

0.1793 






0.69 

0.7238 

0.0870 

0.6593 

0.1606 






0.70 

0.7141 

0.0893 

0.6680 

0.1600 






0.71 

0.7042 

0.0916 

0.6767 

0.1610 






0.72 

0.6940 

0.0939 

0.6854 

0.1421 






0.73 

iO.6835 

0.0963 

0.6939 

0.1346 






0.74 

0.6728 

0.0987 

0.7025 

0.1252 
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Table 4 


The corresponding values of M, B, T, Q/bq, ^ 
for supersonic values of H 


Table 4 (cont’d) 


M 

B 

T 

1/ao 


9/^0 

U 

B 

L 


1.00 

0.0000 

0.1667 

0.9129 

0.0000 

1.00 

1.1180 

0.5000 

0,8000 

0.0350 

1.01 

0.1418 

0.1695 

0.9805 

0.0008 

1.05 

1.1892 

0.6437 

0.8806 

0.0677 

1.02 

0.2010 

0.1722 

0,9280 

0.0080 

1.10 

1.8635 

0.7728 

0.2420 

0.1058 

1.03 

0.8468 

0.1750 

0.9355 

0.0039 

1.15 

1.3409 

0.8934 

0.8645 

0.1483 

1.04 

0.2867 

0.1770 

0.9430 

0.0061 

1.80 

1.4821 

1.0112 

0.2880 

0.1945 

1.05 

0.3202 

0.1007 

0.9504 

0.0085 

1.25 

1.B076 

1.1888 

0.3160 

0.2440 

1.06 

0.3B16 

0.1035 

0.9578 

0.0110 

1.30 

1,5978 

1.2462 

0.3380 

0.2965 

1.07 

0.3807 

0.1863 

0.9652 

0.0140 

1.35 

1.6935 

1.3667 

0.3646 

0.3519 

i.oe 

0.4079 

0.1892 

0.9725 

0.0169 

1.40 

1.7966 

1.4912 

0.3920 

0,4100 

1.09 

0.4337 

0,1980 

0.9798 

0.0200 

1.45 

1.9048 

1.6211 

0.4205 

0.4709 

1.10 

0.4583 

0.1949 

0.9870 

0.0235 

1.50 

2.0226 

1.7581 

0.4500 

0.5345 - 

1.11 

0.4818 

0.1077 

0.9941 

0.0268 

1.B6 

2.1505 

1.8038 

0.4805 

0.5556 

1.12 

0.5044 

0.2006 

1.0014 

0.0304 

1.60 

2.2904 

8.0606 

0.5120 

0.6707 

1.13 

0,5268 

0.8054 

1.0086 

D.0339 

1.65 

8.4448 

8,8309 

0.6445 

0.7436 

1.14 

0.5474 

0.2063 

1.0156 

0.0376 

1.70 

2.6170 

2.4184 

0.5780 

0.8201 

1.15 

0.5679 

0.2092 

1.0226 

0.0415 

1.75 

2.8113 

8.6274 

0.6125 

0.9006 

1.16 

0.5879 

0.2120 

1.0297 

0.0455 

1,00 

5.0339 

2.8644 

0.6480 

0.9866 

1.17 

0.6074 

0.2149 

1.0367 

0.0494 

1.85 

3.2936 

3.1381 

0.6845 

1.0759 

1.18 

0.6264 

0.2178 

1.0436 

0.0535 

1.90 

3.6035 

3.3620 

0.7820 

1.1440 

1.19 

0.6451 

0.2207 

1.0505 

0.0577 

1,96 

3.9346 

3.8570 

0.7605 

1.2768 

1,20 

0,6653 

0.2286 

1.0574 

0.0618 

2.00 

4.4721 

4.3589 

O.BOOO 

1.3900 

1.25 

0,7500 

0.2381 

1.0911 

0.0845 

8.06 

5.1330 

5.0347 

0.8405 

1.5180 

1.30 

0.8307 

0.8586 

1.1288 

Q.1Q7B 

2.10 

6.1133 

6.0310 

0.0880 

1.6634 

1.35 

0.9069 

0.2671 

1.1657 

. 0.1319 

2.15 

7.0247 

7.7605 

0.0245 

1.8401 

1.40 

0,9798 

0.2816 

1.1B66 

0.1567 

2.20 

12.2984 

12.2577 

0.9680 

2.0826 

1.45 

1,0500 

0.2900 

1.2166 

0.1819 

(5 

oo 


1.0000 


1.50 

1.1180 

0.8104 

1.2457 

0.2078 






1.56 

1.1648 

0.3257 

1.8739 

0.2382 






1.60 

1.2490 

0.3336 

1.3018 

0.8590 






1.65 

1.3184 

0-3525 

1.3277 

0.8779 






1.70 

1,3748 

0.3663 

1.3533 

0.3108 






1.75 

1.4361 

0.3798 

1.3781 

0.3376 






1.80 

1.4967 

0.3932 

1.4021 

0.3618 






1.85 

1.5564 

0.4064 

1.4254 

0.3859 
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Table 5 


Auxiliary Functions for Computation of 


H 

M 

— 

F 


\ F^ dA 
‘X'S->‘oO 

dA 

dT 

0.00 

1.0000 

- oQ 

+ 00 

+ CO 

0.0000 

0.05 

0.9998 

-9.5664x10'^ 

1.6676x10^^ 

8.1489x10^^ 

-0.0021 

0.10 

0.9950 

-1.4791x10® 

3.1725x10® 

1.5661x10^° 

-0.0084 

0.15 

0.9887 

-1.2751x10® 

8.0345x10'^ 

3.9617x10® 

-0.0193 

0.20 

0.9798 

-2. 2109x10 4 

5.8077x10® 

2.8547x10 7 

-0.0350 

0.25 

0.9682 

-5.5969x10® 

7.4153x10® 

3.6210x10® 

-0.0561 

0.30 

0.9539 

-1.7922x10® 

1.5497x10® 

6.5149x10® 

-0.0837 

0.55 

0.9368 

-6.7198x10^ 

3.1214x10^ 

1.4797x10® 

-0.1188 

0.40 

0.9155 

-2.8144x10^ 

8.5565x10® 

3,9545x10^ 

-0.1631 

0.45 

0.8931 

-1.2759x102 

2.6555x10® 

1.1903x10^ 

-0.2190 

0.50 

0.8660 

-6.1200x10 

9.O295x102 

3.9424x10® 

-0,2899 

0.65 

0.8352 

-3.0505x10 

3.2815x102 

1.4400x10® 

-0,3806 

0.60 

0.8000 

-1.5558x10 

1.2482x10^ 

6.0566x10^ 

-0.4987 

0.65 

0.7599 

-7.9987 

4.8769x10 

3.1232x10^ 

-0.6559 

0.70 

0,7141 

-4.0789 

1.9195x10 

2. 0784x10 2 

-0.8719 

0.75 

0.6615 

-2.0215 

7.4369 

1.6904x102 

-1.1823 

0.80 

0.6000 

-0.9453 

2.7472 

1.5337x10 2 

-1.6684 

0.85 

0.5267 

-0.3964 

1.1827 

1.4562x10 2 

-2.4667 

0.90 

0.4358 

-0.1336 

0.2477 

1.4039x10 2 

-4.1071 

0.95 

0.3120 

-0.0267 

0.0383 

1.5687x102 

-9.0796 

1.00 

0.0000 

+0.0000 

0.0000 

1.3150x10 2 

ao 
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Table 6 



The. values 

Of df 

d 

and fM 


= 0, 1, 

2, 3, 4. 


df 

d 

f(j^) 

lOf 

fC2] 

fC3] 

fCi! 

0.0000 

20.0000 

1.0000 

0.00000 

0.0000 

0.0000 

0.0000 

0.0025 

20.5000 

1.0506 

0.00256 

0.0003 

0.0000 

0.0000 

0.0050 

21.0000 

1.1025 

0.00525 

0.0012 

0.0000 

0.0000 

0.0100 

22.0000 

1.2100 

0.01103 

0.0053 

0,0002 

0.0000 

0.0150 

23.0000 

1.3225 

0.01736 

0.0124 

0.0006 

0.0000 

0 . 0200 

24.0000 

1.4400 

0.02426 

0.0228 

0.0015 

0.0000 

0.0300 

26,0000 

1.6900 

0.03990 

0 . 0546 

0.0052 

0.0003 

0.0400 

28.0000 

1.9600 

0.05813 

0.1034 

0.0129 

0.0012 

0.0500 

30.0000 

2.2500 

0.07917 

0.1718 

0.0265 

0.0031 

0.0600 

32.0000 

2.5600 

0.10320 

0.2628 

0.0481 

0.0068 

0.0700 

34.0000 

2.8900 

0 . 15043 

0.3793 

0.0799 

0.0131 

0.0800 

36 . 0000 

3.2400 

0.16107 

0.5248 

0.1249 

0.0232 

0 . 0900 

38.0000 

3.6100 

0.19530 

0.7026 

0.1852 

0.0386 

0.1000 

40.0000 

4.0000 

0.23333 

0.9166 

0.2666 

0.0611 

0 . 1100 

42.0000 

4.4100 

0.27537 

0.1170 

0.3706 

0.0927 

0.1200 

44.0000 

4.8400 

0.32160 

0.1468 

0.5022 

1 

0.1361 
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Table 7 

The Values of u = 1 , 2 , 3,4 


A * 

B 

H 

£( 1 ) 

6 ( 2 ) 

£ T ( 3 ) 

£W 

0.0000 

1.0000 

1,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0025 

1 . 0040 

0.9970 

- 0.0007 

- 0.0180 

- 0.1700 

- 2.5000 

0.0050 

1.0080 

0.9940 

- 0.0012 

- 0.0360 

- 0.3300 

- 4.8000 

0.0075 

1.0118 

0 . 9910 

- 0.0016 

- 0.0530 

- 0.4800 

- 6.8000 

0 . 0100 

1.0156 

0.9880 

- 0.0019 

- 0.0700 

- 0.6200 

- 8.6000 

0.0150 

1.0230 

0.9825 

- 0.0027 

- 0.1000 

- 0.8700 

- 11.6500 

0.0200 

1.0300 

0.9757 

- 0.0032 

- 0.1290 

- 1.1300 

- 14.7500 

0.0250 

1.0365 

0.9700 

- 0.0036 

- 0.1560 

- 1.4000 

- 17.6600 

0.0300 

1.0427 

0,9650 

- 0 . 0040 

- 0.1880 

- 1.6700 

- 20.6000 

0.0400 

1.0560 

0.9545 

- 0.0047 

- 0.2400 

- 2.1100 

- 25.2000 

0.0500 

1.0690 

0 . 9440 

- 0.0042 

- 0.2910 

- 2.4700 

- 29.2700 

0.0600 

1.0813 

0.9325 

- 0.0037 

- 0,3300 

- 2.7900 

- 33.1000 

0.0700 . 

1.0945 

0.9210 

- 0.0023 

- 0.3710 

- 3.1500 

- 36.4000 

0.0800 

1.1075 

0.9115 

- 0.0010 

- 0.4090 

- 3.4500 

- 39.4500 

0.0900 

1.1210 

0.9010 

+ 0 . 0007 

- 0.4430 

- 3.7000 

- 42.1600 

0.1000 

1.1335 

0.8915 

0.0028 

- 0.4740 

- 3,9400 

- 44.1900 

0 . 1100 

1.1467 

0.8825 

0.0053 

- 0.5060 

- 4.1200 

- 46.1000 

0.1200 

1 . 1610 

0.8735 

0.0080 

- 0.5310 

- 4.3000 

- 48.0000 

0.1300 

1.1755 

0.8642 

0.0110 

- 0.5570 

- 4.4630 

- 49.5000 

0 . 1400 

1,1850 

0.8554 

0.0143 

- 0.5840 

- 4.6150 

- 50.8000 

0.1500 

1.1972 

0.8476 

0.0181 

- 0.6100 

- 4.7500 

- 52.0300 

0.1600 

1.2110 

0.8390 

0.0219 

- 0.6310 

- 4.9000 

- 53.0500 

0.1700 

1.2235 

0.8295 

0,0260 

- 0.6500 

- 5.0800 

- 54.0000 

0.1800 

1.2368 

0.8208 

0.0300 

- 0.6700 

- 5.1400 

- 54.9000 

0.1900 

1 . 2496 

0.8120 

0.0349 

- 0.6880 

- 5.2450 

- 55.7200 

0.2000 

1.2615 

0.8040 

0.0398 

- 0.7060 

- 5.3370 

- 56.6200 

0 . 2100 

1,2732 

0.7960 

0.0443 

- 0.7250 

- 5.4260 

- 57.4000 

0 . 2200 

1.2856 

0.7885 

0.0493 

- 0.7420 

- 6.5100 

- 58.0800 

0 . 2300 

1.2980 

0.7810 

0.0547 

- 0.7560 

- 5.5950 

- 58.5500 

0.2400 

1.3117 

0.7740 

0.0603 

- 0.7710 

- 5,6720 

- 59.0000 

0.2500 

1.3255 

0.7664 

0.0668 

- 0.7860 

- 5.7400 

- 59.4500 

0.2600 

1.3370 

0.7585 

0.0727 

- 0.8000 

- 5.8050 

- 59.8500 
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0.0025 

0.0000 

0.0025 

0.00125 

0.3765 

1 

40.3770 

+153.2134 

0.1915 


0.0025 

0.0050 

0.00000 

0.0000 

40.7513 


0.0000 

0.0050 

0.0000 

0.0025 

0.0050 

0.0075 

0. Q 0250 

0.00125 

0.7488 

0.3722 


m-m 



0.0050 

0.0100 

0.00000 

0.0000 

41.4911 

161.1307 

0.0000 

0.0100 

0.0000 

0.0050 

! 0.0100 
[ 0.0150 

0.005 QD 

0.00250 

1.4801 

0.7317 

41.4910 

42.2397 


0.8046 

0.4160 


0.0100 

0.0200 

0.00000 

0.0000 

42.9203 

171.4925 

0.0000 

0.0150 

0.0000 

0.0150 

0.00750 

2.1950 

42.2395 

165.9949 

1.2449 


0.0050 

0.0200 

0.00500 

1.4447 

42.9202 

171.2956 

0.8564 


O.OIOQ 

0.0250 

0.00250 

0.7138 

43.6366 

176.5320 

0.4413 


0.0150 

0.0300 

0.00000 

0.0000 

44.3730 

181.7387 

0.0000 

0.0200 

0.0000 

0.0200 

0.01000 

2.8893 

42.9198 

170.7048 

1.7071 


0.0050 

0.0250 

0.00750 

2.1415 

43.6363 

176.1484 

1.3211 


0.0100 

0.0300 

0.00500 

1.4116 

44.3728 

181.5518 

0.9077 


0.0200 

0.0400 

0.00000 

0.0000 

45.7909 

192.4503 

0.0000 

0.0300 

0.0000 

0.0300 

0.01500 

4.2344 

44.3713 

180.0566 

2.7010 


0.0100 

0.0400 

0.01000 

2.7592 

45.7899 

191.7377 

1.9174 


0.0200 

0.0500 

0.00500 

1.3483 

47.1679 

203.1159 

1.0155 


0.0300 

0.0600 

0.00000 

0.0000 

48.4506 

213.7738 

0.0000 

0.0400 

0.0100 

0.0500 

0.01500 

4.0446 

47.1655 

201.7565 

3.0264 


0.0200 

0.0600 

0.01000 

2.6317 

48.4491 

213.1260 

2.1313 


0.0300 

0.0700 

0.00500 

1.2846 

49.6875 

224.6336 

1.1231 

j 

0.0400 

0.0800 

0.00000 

0.0000 

50.9896 

236.4134 

0.0000 

0.0600 

0.0300 

0.0900 

0.01500 

3.6830 

52.1909 

246.9304 

3.7041 


0.0400 

0.1000 

0.01000 

2.4008 

53.4178 

259.4997 

2.5951 

1 

1 

0.0500 

0.1100 

0.00500 

1.1746 

54.6372 

272.4613 

1.3623 


0.0600 

0.1200 

0.00000 

0.0000 

55.8202 

285.9281 

0.0000 

0.0800 

0.0500 

0.1300 

0.01500 

3.3715 

56.9404 

297.3904 

4.4610 


0.0600 

0.1400 

0.01000 

2.1978 

58.0675 

310.1456 

3.1016 


0.0700 

0.1500 

0.00500 

1.0759 

59.2298 

324.0247 

1.6201 


0.0800 

0.1600 

0.00000 

0.0000 

60.3022 

338.3378 

0.0000 

0.1000 

0.0800 

0.1800 

0.01000 

2.0111 

62.2637 

365.5425 

3.6556 


0.0900 

0.1900 

0.00500 

0.9830 

63.2227 1 

380.0431 

1.9002 


0.1000 

0.2000 

0.00000 

0.0000 

64.2081 

394.3472 

0.0000 

0.1100 

0.0900 

0.2000 

0.01000 

1.9229 

64.2037 

393.9469 

3.9396 


0.1000 

0.2100 

0.00500 

0.9401 

65.1558 

408.4695 

2.0423 


0.1100 

0.2200 

0.00000 

0.0000 

66.1189 

423.7961 

0.0000 

0.1200 

0.1000 

0.2200 

0.01000 

1.8391 

66.1141 

423.4159 

4.2343 


0.1100 

0.2300 

0.00500 

0.8996 

67.0534 

439.1157 

2.1956 


0.1200 

0.2400 

0.00000 

0.0000 

68.0051 

456.3483 

0.0000 
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The Values of M ^, 

“ 2 ^ “ 3 » 

“4 


r 


2 A * = 

m 

“2 

“3 

“4 

0.0025 

0.0000 

0.0025 

20.4992 

- 0.0005 

0.0000 

- 19.9993 


0.0025 

0.0050 

20.4986 

- 0.0005 

0.0005 

- 20.4986 

0,0050 

0.0000 

0.0050 

20.9984 

- 0.0011 

0.0000 

- 19.9988 


0.0025 

0.0075 

20.9979 

- 0.0010 

0.0004 

- 20.4981 


0.0050 

0.0100 

20.9975 

- 0.0009 

0.0009 

- 20.9975 

0.0100 

0.0000 

0.0100 

21 . 9968 

- 0.0021 

0.0000 

- 19.9981 


0.0050 

0.0150 

21.9955 

- 0.0018 

0.0007 

- 20.9964 


0.0100 

0.0200 

21.9946 

- 0.0014 

0.0014 

- 21.9946 

0.0150 

0.0000 

0.0150 

22.9945 

- 0.0031 

0.0000 

- 19.9973 


0.0050 

0. 0200 

22.9933 

- 0.0025 

0.0006 

- 20.9957 


0.0100 

0.0250 

22.9923 

- 0.0020 

0 . 0011 

- 21.9938 


0.0150 

0.0300 

22.9912 

- 0.0015 

0.0015 

- 22.9912 

0.0200 

0.0000 

0.0200 

23.9919 

- 0.0039 

0.0000 

- 19.9968 


0,0050 

0.0250 

23.9906 

- 0.0031 

0.0004 

- 20.9951 


0.0100 

0.0300 

23.9892 

- 0.0024 

0.0008 

- 21.9929 


0.0200 

0.0400 

25.9868 

- 0.0011 

0.0011 

- 25.9868 

0.0300 

0.0000 

0.0300 

25.9847 

- 0.0051 

0.0000 

- 19.9960 


0.0100 

0.0400 

25.9811 

- 0.0029 

0.0002 

- 21.9915 


0.0200 

0.0500 

26.9797 

- 0.0008 

- 0,0000 

- 23.9862 


0.0300 

0.0600 

25.9788 

+ 0.0009 

- 0.0009 

- 25.9788 

0.0400 

0.0100 

0 . 0500 

27.9719 

- 0.0025 

- 0.0002 

- 21.9915 


0.0200 

0.0600 

27.9702 

+ 0.0001 

- 0.0010 

- 25.9859 


0.0300 

0.0700 

27.9701 

+ 0.0030 

- 0,0028 

- 25.9793 


0.0400 

0,0800 

27.9701 

+ 0.0050 

- 0.0056 

- 27.9701 

0.0600 

0.0300 

0.0900 

31.9443 

+ 0.0108 

- 0.0061 

- 25.9812 


0.0400 

0.1000 

31.9457 

+ 0.0147 

- 0.0100 

- 27.9732 


0.0500 

0.1100 

31.9485 

+ 0.0185 

- 0.0153 

- 29.9635 


0.0600 

0.1200 

31.9520 

+ 0.0215 

- 0.0215 

- 31.9520 

0.0800 

0.0500 

0.1300 

35.9163 

+ 0.0356 

- 0.0204 

- 29.9716 


0.0600 

0,1400 

36.9217 

+ 0.0405 

- 0.0284 

- 31.9617 


0.0700 

0.1500 

35.9289 

+ 0.0454 

- 0.0382 

- 33.9505 


0.0800 

0.1600 

35 . 9369 

+ 0.0499 

- 0.0499 

- 35.9369 

0.1000 

0.0800 

O.ISOO 

39.9019 

0.0794 

- 0.0678 

- 35.9554 


0.0900 

0.1900 

39.9161 

0.0853 

- 0.0731 

- 37.9444 


0.1000 

0.2000 

59.9304 

0.0907 

- 0.0907 

- 39.9304 

0.1100 

0.0900 

0.2000 

41.8976 

0.1048 

- 0.0775 

- 37.9578 


0.1000 

0 . 2100 

41.9109 

0.1114 

- 0.0962 

- 39.9430 


0.1100 

0.2200 

41.9270 

0.1169 

- 0.1169 

- 41.9270 

0.1200 

0.1000 

0.2200 

43.8898 

0.1342 

- 0.1007 

- 39.9580 


0.1100 

0.2300 

43.9100 

0.1407 

- 0.1223 

- 41.9458 


0.1200 

0.2400 

43.9309 

0.1472 

- 0.1472 

- 43.9309 
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The Values of V ^, Vg , V , 


f* 

f 

2 A * = 

J*M 

* 

B 


V 

f* 

0.0025 

0.0000 

0.0025 

1 

.0506 

- 1.0000 

0.0506 

0.0504 


0.0025 

0.0050 

0 

.0000 

0.0000 

0.0000 

0.0000 

0.0050 

0.0000 

0.0050 

1 

.1024 

- 0.0000 

0.1024 

0.1018 


0.0025 

0.0075 

1 

.1024 

- 1.0506 

0.0518 

0.0514 


0.0050 

0.0100 

0 

.0000 

0.0000 

0.0000 

0.0000 

0.0100 

0.0000 

0.0100 

1 

.2099 

- 1.0000 

0.2099 

0.2074 


0.0050 

0.0150 

1 

.2099 

- 1.1024 

0.1074 

0.1056 


0.0100 

0.0200 

0 

.0000 

0.0000 

0.0000 

0.0000 

0.0150 

0.0000 

0.0150 

1 

.3224 

- 1.0000 

0.3224 

0.3168 


0.0050 

0.0200 

1 

.3224 

- 1.1024 

0.2199 

0.2146 


0.0100 

0.0250 

1 

.3224 

- 1.2099 

0.1124 

0.1091 


0.0150 

0.0300 

0 

.0000 

0.0000 

0.0000 

0.0000 

0.0800 

0.0000 

0.0200 

1 

.4398 

- 1.0000 

0.4398 

0.4292 


0.0050 

0.0250 

1 

.4398 

- 1.1024 

0.3373 

0.3272 


0.0100 

0.0300 

1 

.4398 

- 1.2099 

0.2299 

0.2218 


0.0800 

0.0400 

0 

.0000 

0.0000 

0.0000 

0.0000 

O.OSOO 

0.0000 

O.OSOO 

1 

.6897 

- 1.0000 

0.6897 

0.6655 


0.0100 

0.0400 

1 

.6896 

- 1.2099 

0.4797 

0.4579 


0.0200 

0.0500 

1 

.6896 

- 1.4398 

0.2498 

0.2358 


0.0300 

0.0600 

0 

.0000 

0.0000 

0.0000 

0.0000 

0.0400 

0.0100 

0.0500 

1 

.9594 

- 1.2099 

0.7494 

0.7075 


0.0200 

0.0600 

1 

.9594 

- 1.4398 

0.5195 

0.4845 


0.0300 

0.0700 

1 

.9594 

- 1.6896 

0.2697 

0.2484 


0.0400 

0.0800 

0 

.0000 

0.0000 

0.0000 

0.0000 

0 . 0600 

0.0300 

0.0900 

2 

.5687 

- 1,6897 

0.8689 

0.7829 


0.0400 

0.1000 

2 

t.5588 

- 1.9596 

0.5992 

0.5341 


0.0500 

0.1100 

2 

.5589 

- 2.2494 

0.3095 

0.2731 


0.0600 

0.1200 

C 

.0000 

0.0000 

0.0000 

0.0000 

0.0800 

0.0500 

0.1300 

*2 

>.2381 

- 2.2497 

0.9883 

0.8541 


0.0600 

0.1400 

*2 

c. 

>.2385 

- 2.5596 

0.6788 

0.5806 


0.0700 

0.1500 

•2 

>.2390 

- 2.8895 

0.3494 

0,2961 


0.0800 

0.1600 

C 

► .0000 

0.0000 

0.0000 

0,0000 

0.1000 

0.0800 

0.1800 

*5 

5.9991 

- 3.2405 

0.7586 

0.6226 


0.0900 

0.1900 

4.0001 

- 5.6107 

0.3893 

0.3161 


0.1000 

0.2000 

C 

).0000 

0.0000 

0.0000 

0.0000 

0.1100 

0.0900 

0.2000 

4 

L.4101 

- 5.6116 

0.7985 

0.6420 


0.1000 

0.2100 

4.4111 

- 4.0018 

0.4092 

0.3257 


0.1100 

0.2200 

C 

5.0000 

0.0000 

0.0000 

0.0000 

0.1200 

0.1000 

0.2200 

4.8415 

- 4.0028 

0.8385 

0.6611 


0.1100 

0.2300 

4 

t.8428 

- 4.4135 

0.4292 

0.3352 


0.1200 

0.2400 

0.0000 

0.0000 

0.0000 

0.0000 


83 









NAOA TN Ho. 1096 
Table 11 

Streamlines of the Supersonic Flow (53) in the Physical Plane 


Note that x and y are given only up to a multiplicative 
constant and that the flow is symmetric with respect to the 
x-axis. 



0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


II 

= 0,1 

X 

y 

4.18 

0.45 

8.54 

0.46 

13.78 

0.47 

19,58 

0.48 

26.70 

0.49 

34.70 

0.50 



4.03 

8.19 

13.21 

19.03 

25.91 

33.81 


0.90 

0.92 

0.94 

0.96 

0.98 

0.99 



4,11 

8.39 

13.39 

19.29 

26.09 

33.79 



1.46 

1.49 

1.53 


4.00 

8.12 

13.12 

18.74 

25.68 

33.18 


1.82 

1.86 

1.90 

1.94 

1.98 

2.03 



4.06 

8.32 

13.26 

20.16 

26.40 

34.60 


2.27 

2.32 

2.37 

2.42 

2.48 

2.54 


4.28 

8.42 

15.42 

19.16 

25.96 

34.12 


2.72 

2.78 

2.86 

2.39 

2.96 

3.02 
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